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Abstract: We further develop the formalism for taking position-space cuts of eikonal 
diagrams introduced in ref. [1]. These cuts are applied directly to the position-space rep¬ 
resentation of any such diagram and compute its discontinuity to the leading order in the 
dimensional regulator. We provide algorithms for computing the position-space cuts and 
apply them to several two- and three-loop eikonal diagrams, finding agreement with results 
previously obtained in the literature. We discuss a non-trivial interplay between the cutting 
prescription and non-Abelian exponentiation. We furthermore discuss the relation of the 
imaginary part of the cusp anomalous dimension to the static interquark potential. 
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1 Introduction 


The infrared singularities of gauge theory scattering amplitudes play a fundamental role 
in particle physics for phenomenological as well as more theoretical studies. Determining 
the long-distance singularities is necessary for combining the real and virtual contributions 
to the cross section, as the divergences of the separate contributions only cancel once they 
are added. Infrared singularities moreover dictate the structure of large logarithmic contri¬ 
butions to the cross section, allowing such terms to be resummed—which is in many cases 
required in order to obtain reliable perturbative predictions. Beyond their significance to 
collider phenomenology, long-distance singularities are highly interesting from a theoretical 
point of view. Among several properties, they have a universal structure among different 
gauge theories; moreover, their exponentiation properties [2-14] and their relation to the 
renormalization of Wilson line correlators [15-21] allow their perturbative expansion to be 
explored to all orders, a feat currently unattainable for complete scattering amplitudes. 

The basic tool for computing the infrared singularities of any scattering amplitude is 
provided by the eikonal approximation. In this limit the momenta of the soft gauge bosons 
emitted between the partons emerging from the hard interaction are neglected with respect 
to the hard momenta pi . As a result, each hard parton i simply acts as a source of soft 
gluon radiation and is accordingly replaced by a semi-infinite Wilson line 


( I\;.,; = V exp 


(ig J d tVi- , 


( 1 . 1 ) 


which extends from time t = 0 , when the hard scattering takes place, to infinity along the 
classical trajectory of the hard parton, traced out by its four-velocity . The long-distance 
singularities of the scattering amplitude of the hard partons are then encoded in the eikonal 
amplitude 

S{p/ij i e) = (01$.^ (8) <g> $<J0), (1.2) 

which has the same soft singularities as the original amplitude, but is much simpler to 
compute. An important feature of the eikonal amplitude (1.2) is the fact that it depends 
on the kinematics only through the angles 7 ij between the four-velocities (defined through 
cosh 7 ^ = | Vi ■ Vj\). Before renormalization, the integrals involved in the loop-level contri¬ 
butions to S are thus scale invariant and vanish identically. This in turn allows the infrared 
singularities at any loop order to be computed by studying the ultraviolet renormalization 
factor of the Wilson line correlator (1.2) [19, 21-25]. This renormalization factor forms a 
matrix in the space of color configurations available for the scattering process at hand, 
referred to as the soft anomalous dimension matrix. In processes involving only two Wilson 
lines, this matrix reduces to the cusp anomalous dimension, a quantity which has been 
computed in QCD up to three loops [21, 26, 27]. In M = 4 super Yang-Mills theory, the 
cusp anomalous dimension is known to three loops [28], and partial results have been ob¬ 
tained at four loops [29, 30]. For multi-parton amplitudes, the soft anomalous dimension 
matrix has been computed through two loops for massless [31, 32] as well as massive [33- 
37] Wilson lines. Recently, much progress has been made toward the calculation of the soft 
anomalous dimension matrix at three loops [38, 39]. 
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In this paper we continue exploring a notion of cuts of eikonal diagrams (i.e., the dia¬ 
grams contributing to the eikonal amplitude) introduced in ref. [1]. Applied to any eikonal 
diagram, the cuts compute the discontinuities of the diagram, in analogy with the Cutkosky 
rules for standard Feynman diagrams. The discontinuities are in turn readily combined 
to produce the imaginary part of the diagram, a direct computational method of which 
is desirable in several contexts. Indeed, collinear factorization theorems for non-inclusive 
observables were pointed out in refs. [40, 41] to be violated due to exchanges of Glauber- 
region (i.e., maximally transverse) gluons. The resulting factorization-breaking terms are 
purely imaginary and take the form of the non-Abelian analog of the QED Coulomb phase. 
Therefore, by utilizing the all-order exponentiation property of the eikonal amplitude, the 
latter could be obtained directly by computing the imaginary part of the exponent. The 
resulting non-Abelian Coulomb phase [42, 43] may also aid studies of interference effects. 
The importance of understanding the imaginary part of eikonal diagrams has also recently 
been highlighted in studies regarding rapidity gaps [44, 45] . Moreover, cuts of Wilson line 
correlators are naturally relevant for cross section calculations [46, 47]. 

A cutting prescription for eikonal diagrams may also provide the first step toward ex¬ 
tending the modern unitarity method [48-57] to eikonal amplitudes. The development of 
the unitarity method has led to a dramatic improvement in the ability to compute loop- 
level (non-eikonal) scattering amplitudes at high multiplicity. In this approach, the loop 
amplitude is decomposed into a linear basis of loop integrals which are computed inde¬ 
pendently (for example, by means of Feynman parametrization, or differential equations 
[58, 59]). The calculation of the loop amplitude is then reduced to the problem of determin¬ 
ing the integral coefficients. This step is performed by applying to both sides of the basis 
decomposition of the loop amplitude a number of cuts which have the effect of putting the 
internal lines on shell. In basic unitarity (as opposed to generalized unitarity), the cuts 
employed measure the discontinuity of the amplitude in its various kinematical channels. 
Unitarity has proven highly successful, notably in computing one-loop amplitudes with 
many partons in the final state. It is therefore natural to look for extensions of this method 
to other physical quantities with a perturbative expansion. 

It should be emphasized that Cutkosky rules for eikonal diagrams have been introduced 
previously in the literature, as a cut prescription applied directly to the momentum-space 
representation of the diagrams [21]. In contrast, the cuts introduced in ref. [1] and further 
studied here are applied to the position-space representation of the eikonal diagrams. A 
notion of position-space cuts of non-eikonal diagrams exists in the literature in the form of 
a cutting equation that follows from Veltman’s largest-time equation [60]. However, that 
notion is conceptually different from the position-space cuts in this paper, since the former 
has the effect of cutting a given diagram in two disconnected subdiagrams while the latter 
does not. Moreover, in practice, the largest-time equation is typically not applied directly, 
but rather serves to derive the momentum-space Cutkosky rules, which in turn are used 
to obtain the imaginary part of a diagram. As already observed in ref. [1], position-space 
cuts provide a substantial simplification over momentum-space cuts in the computation of 
imaginary parts of eikonal diagrams. There has been recent interest in the literature in 
studying Wilson line correlators in position space, in particular refs. [61-63] which inves- 
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tigate the structure of infrared singularities and factorization in position space. Moreover, 
position-space analogs of generalized unitarity cuts of Wilson line correlators were recently 
introduced in ref. [64]. 

The structure of this paper is as follows. In section 2 we discuss the origin of the 
imaginary part of Wilson line correlators from the point of view of causality as well as 
unitarity. We then show how the imaginary part can be computed from the position- and the 
momentum-space representations at one loop. In section 3 we review the formula in ref. [1] 
for the imaginary part of L-loop eikonal diagrams containing no internal (i.e., three- or four- 
gluon) vertices to the leading order in the dimensional regulator e. We furthermore discuss 
the relation of the discontinuities of the diagrams to their imaginary part. In section 4 
we apply the formalism to compute the imaginary part of a number of two- and three- 
loop diagrams and discuss a non-trivial interplay between the cutting prescription and 
non-Abelian exponentiation. In section 5 we turn to formulas for the imaginary part of 
eikonal diagrams with internal vertices and provide details on its computation. We give 
our conclusions in section 6. Appendix A explains our method for computing the principal- 
value integrals involved in the cutting prescription. In appendix B we present our algorithm 
for re-expressing multiple polylogarithms in terms of ones with constant indices. 

2 Imaginary parts of eikonal diagrams and their physical origin 

In this section we will discuss the origin of the imaginary part of Wilson line correlators from 
the point of view of causality as well as unitarity. These viewpoints are naturally provided 
by the representation of the correlator in position and momentum space, respectively. We 
will show how the imaginary part can be computed directly from each of the respective 
integral representations at one loop. 

We adopt the convention that all velocities are outgoing, such that the velocities as¬ 
sociated with outgoing and incoming states respectively have positive and negative time 
components. We will take the gauge group to be SU(N) and work in Feynman-’t Hooft 
gauge with (+, —, —, —) spacetime signature. Ultraviolet divergences will be regulated by 
computing all diagrams in D = 4 — 2e dimensions with e > 0. To avoid complications 
arising from regulating collinear singularities, we take all velocities to be time-like, v'f = 1 . 


We start our investigations by examining the simplest eikonal diagram, the one-loop 
exchange illustrated in figure 1. In both kinematic regions (a) and (b), the position-space 
representation of the diagram is straightforwardly obtained by direct perturbative expan¬ 
sion 1 in g of the correlator (1.2) and takes the form 

f ( i) =d r du r dt 2 ——, (2.D 

Jo Jo [~(hvi — t 2 v 2 ) 2 + irj\ 

where the prefactor is defined as C ^ = g 2 C , jr r ^'^ /2 1 ^ , with Cf = N 2 n X denoting the 
quadratic Casimir of the fundamental representation. Furthermore, t\.t 2 have the dimen- 

1 See chapter 8 of ref. [65] for the position-space Feynman rules. 
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(a): V\ ■ V 2 > 0 (b): v\ • V 2 < 0 

Figure 1. One-loop eikonal diagrams. In (a) both Wilson lines represent final-state partons (as for 
example in e + e _ —> qq ). Thus, the external velocities are in the region v\ ■ V 2 > 0, and the partons 
can become lightlike separated. In (b) one Wilson line represents a final-state parton, and one 
represents an initial-state parton (as for example in deep inelastic scattering). Thus, the external 
velocities are in the region V\ ■ i) 2 < 0, and the partons are never lightlike separated. 


sion of time and denote the positions of the attachment points of the soft-gluon propagator 
on the Wilson lines spanned by the four-velocities v± and V 2 - 

The integrations in eq. (2.1) produce an infrared divergence which can be extracted via 
the change of variables (ti, ^ 2 ) = (Ax, A(1 — x)) with 0 < x < 1, where A has the dimension 
of length, 





dx- 


V\ ■ V2 


'0 [— (xvi — (1 — x)v 2 ) 2 + il]\ 


l—e 


( 2 . 2 ) 


Indeed, the A-integral is has an infrared divergence, owing to the exchange of gluons of 
increasingly longer wavelength as A —> 00 . This divergence can be regularized in a gauge 
invariant fashion by introducing an exponential damping factor e -AA with A< 1 , whereby 
it becomes 




2e 


dA e 


-AA 


A 1-2e 


= r(2e) 


= J_ /MA 


A ) 


2e 


2 e VA ) 


+ O(e 0 


(2.3) 


The two diagrams in figure 1 have the same integrand; however, as the external kine¬ 
matics is taken from the distinct regions v\ ■ V 2 > 0 and rq • V 2 < 0 , the integrations will 
produce distinct results. It is most convenient to compute the diagram in figure 1(b) first 
and obtain the result for figure 1(a) by analytic continuation as follows. For the diagram in 
figure 1 (b), we may define the deflection angle 7 > 0 such that cosh 7 = —v\ ■ V 2 , in terms 
of which the diagram in figure 1 (b) becomes, to the leading order in e, 


,(D = C« 

i(b) 2e 



2e 

7 coth 7 . 


(2.4) 


Likewise, for the diagram in figure 1(a), we may define the cusp angle 7 > 0 such that 
cosh 7 = vi ■ V 2 - The integrated expression for this diagram can thus be obtained from 
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eq. (2.4) by replacing 7 ^ -> ni — 7 ^ = 7 ^, 

F 1% = %- e - (f) 2£ ^ - 7r *) coth ^ • ( 2 - 5 ) 

We observe that the imaginary parts of the eikonal diagrams in figs. 1(a) and 1(b) are re¬ 
spectively non-vanishing and vanishing. Before turning to the question of how the imaginary 
parts of the diagrams in eqs. (2.4)-(2.5) may be extracted from their integral representation 
in eq. ( 2 . 2 ), let us consider their physical origin and interpretation. 

From the position-space representation (2.1) of the eikonal diagram, the origin of the 
imaginary part may be understood from a simple causality consideration as follows. As our 
focus is on computing the imaginary part to the leading order in e, the e in the propagator 
exponent may be dropped once the infrared divergence has been extracted. After moreover 
stripping off real prefactors from eq. ( 2 . 1 ), the integral takes the form, 

poo poo 1 

/ dil / dt 2 7,-7---• ( 2 -6) 

Jo Jo {tivi - t 2 v 2 y - iv 

Now, for the kinematics corresponding to the diagram in figure 1(a), there are regions 
= e ±7 within the integration domain where (t\Vi — t 2 v 2 ) 2 = 0. Here the —if] term 
becomes relevant and generates an imaginary part. What is happening physically at such 
times t\,t 2 is that the two partons traveling along v\ and v 2 become lightlike separated. 
This is illustrated in figure 2(a). As a result, the phases of their states will change through 
exchanges of lightlike gluons (or photons)—leading to observable consequences that will 
be discussed shortly. In contrast, for external kinematics corresponding to the diagram 
in figure 1 (b), the integral in eq. ( 2 . 6 ) has a vanishing imaginary part: the denominator 
(tivi —t 2 v 2 ) 2 is strictly positive within the region of integration, and the —iv can therefore 
be dropped. In this situation, the partons are never lightlike separated, as illustrated in 
figure 2 (b), and the phases of their states cannot change through exchanges of lightlike 
massless gauge bosons. 

These observations on the evolution of the phases of the hard-parton states suggest that 
the imaginary part of the correlator of two Wilson lines defines an interparton potential. 
Indeed, in the non-relativistic limit, the final and initial two-particle states are related in 
the interaction picture through time evolution by \ f)j = e l Jo dte AtVlt \i)j where Vj denotes 
the interaction potential. The relation of the correlator to a non-relativistic potential can 
be made precise in the situation where the pair of energetic particles carry no color charges, 
as for example in the case of an e + e _ pair. In Abelian gauge theories, the correlator of two 
Wilson lines can be written as the exponential of the sum of connected diagrams [2], 

W = ($ V1 $ V2 ) = exp (F« + 0(g 4 )) , (2.7) 

where F^ is the one-loop diagram in figure 1 , and the additional diagrams contain a 
single lepton loop connected to the Wilson lines by an arbitrary (even) number of soft- 
photon exchanges. Using the result for the diagram F^ computed for time-like kinematics 
in eq. (2.5) (with Cf —> 1 to recover the Abelian case), the anomalous dimension of the 
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Figure 2. The one-loop eikonal diagrams of figure 1 embedded in a space-time diagram. In (a) both 
Wilson lines are confined to the interior of the future light cone, describing two final state partons. 
In this case the partons can become lightlike separated, which is illustrated by the exchange of a 
lightlike gluon (i.e. a gluon that is aligned to the light cone). In (b) there is one incoming parton 
inside the past light cone and an outgoing parton inside the future light cone. These partons are 
never lightlike separated. Indeed, the gluon stretching between these two Wilson lines is necessarily 
time-like (i.e. off shell). 


Wilson-line correlator—i.e., the cusp anomalous dimension—evaluates to 

r cu , p ( 7 ) ^ - lim = -^(7-ri)»th 7 . (2.8) 

The non-relativistic limit corresponds to the small-angle regime 7 ~ 0 where the two 
velocities v\ and V2 are nearly collinear, and the relative velocity of the hard leptons thus 
small. Accordingly, expanding eq. ( 2 . 8 ) around 7 = 0 and taking the imaginary part, we 
find 

Im r cusp ( 7 ) = £- + 0(7°). ( 2 . 9 ) 

We observe that the imaginary part of the cusp anomalous dimension evaluated in time¬ 
like kinematics takes the form of the non-relativistic Coulomb potential (the appropriate 
dimension of energy is acquired after replacing the angle 7 by the distance between the 
two fermions). 

This relation does not extend to generic non-Abelian gauge theories, as we will discuss 
shortly. It does, however, extend to the case of conformal field theories, such as J\T = 4 super 
Yang-Mills theory, where the state-operator correspondence relates Wilson-line operators 
in Minkowski space to states in M x AdS 3 . In radial quantization, a pair of Wilson lines 
intersecting at a cusp angle 7 with the resulting anomalous dimension T cusp (7) is mapped to 
a pair of static charges in AdS 3 separated by a distance of 7 with an electrostatic energy 2 

2 The real part of the cusp anomalous dimension gives rise to an imaginary part of the electrostatic 
energy. As argued in ref. [66], the resulting non-unitary time evolution is accounted for by the real radiation 
of soft and collinear gluons along the Wilson lines. 
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of Imr cusp [ 66 ]. For small values of the cusp angle, the charges on AdS 3 become closer 
than the curvature scale, and the electrostatic energy takes the form of the non-relativistic 
interquark potential in flat space [28, 67]. (The non-relativistic approximation becomes 
relevant here, as in the small-angle regime 7 ~ 0 , the relative velocity of the hard partons 
is small, as discussed above.) 

However, for non-Abelian and non-conformal gauge theories such as QCD, diagrams 

containing loop corrections to the soft propagators will have a dependence on the beta 

function, thereby explicitly breaking the scale invariance of the diagram. As a result, in 

( 3 ) 

QCD, the imaginary part of the three-loop cusp anomalous dimension Tcusp differs from 
the static interquark potential by terms proportional to the beta function [27]. (This can 
be seen by comparing the Nj contribution to Tcusp, given in eq. (A.2) of ref. [ 68 ], against 
the Nj term of the three-loop static QCD potential 3 , given in eq. (10) of ref. [69].) 


Let us now turn to the question of how the imaginary part of the eikonal diagrams in 
figure 1 may be obtained from their integral representation in eq. ( 2 . 2 ) where the infrared 
divergence has been extracted. We will restrict attention to the leading order in the dimen¬ 
sional regulator e, and accordingly drop the e in the propagator exponent. We can then 
utilize the formula 

L ix WjTT n = pv / dx M T "l a* Hx)6(D(x)) . (2.10) 

where PV indicates that the Cauchy principal value prescription is to be applied, and 
the integration bounds a and b are real numbers. The denominator D(x ) is a real-valued 
polynomial in x, and the numerator f(x) is an arbitrary real-valued function with no poles 
or branch points inside the integration path. As both integrals on the right-hand side of 
eq. ( 2 . 10 ) are real, this formula achieves a decomposition into a purely real and purely 
imaginary part. 

Accordingly, at one loop, we define the position-space cut prescription 


1 

D(x) ± IT ] 


cut 


T? rid(D(x)) , 


( 2 . 11 ) 


in terms of which it is straightforward to obtain the imaginary part of the diagrams in 
figure 1 to the leading order in e. For example, considering the time-like kinematics situation 
in figure 1 (a) and applying the prescription ( 2 . 11 ) to eq. ( 2 . 2 ) with the e in the propagator 
exponent set to zero, we find 

poo j \ p — A \ pi 

ImF^) = -ttC w (vi • v 2 )p 2e J —y[Z 2 —J dx5((xvi - (1 - x)v 2 ) ) • (2.12) 

We can integrate out the delta function by use of the identity 
pi 1 

/ d x 5(Ax 2 + Bx + C) = —j= 9(pi)9(l — pi), (2.13) 

Jo VA i=12 

li Note that in the literature on the interquark potential, the loop order is often defined as one less than 
the standard notion. 
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where A = B 2 — 4AC and pi respectively denote the discriminant and roots of the poly¬ 
nomial. The roots pi = of the delta function argument in eq. (2.12) are manifestly 

located inside the domain of the x-integration. The result of integrating out the delta 
function in eq. (2.12) is therefore 

ImF l ( (a) = (f) C ° th7 ’ (2 ' 14) 

in agreement with eq. (2.5). The calculation for the diagram in figure 1(b) is completely 
analogous, except that in this case dx S((xvi — (l — x)v 2 ) ) = 0, as both roots pi = 1 _^ ±7 
are located outside the domain of integration. We therefore find a vanishing imaginary 
part, in agreement with eq. (2.4). We conclude that in both cases (a) and (b), the cutting 
prescription (2.11) produces the correct imaginary part. We introduce a graphical notation 
for the cutting prescription (2.11) in figure 3. 



Figure 3. Graphical representation of a cut of an eikonal diagram in its position-space representa¬ 
tion. The black dots represent the emission and absorption of a lightlike gauge boson. 


It is natural to ask whether the imaginary part of eikonal diagrams can also be obtained 
from their momentum-space representation 1 


F w = ig 2 C F p 2e 


d D k v\ ■ V 2 

(2 n) D ( k 2 + ivj){yi ■ k + irj)(v2 ■ k — irj) 


(2.15) 


Such a cutting prescription was provided in ref. [21]. Here it was shown that the imaginary 
part of the one-loop diagram in eq. (2.15) may be obtained by replacing the two eikonal 
propagators by delta functions, 


1 

k ■ Vi± ip 


cut 


=p27r i6(v®)5(k ■ Vi). 


(2.16) 


This prescription can be thought of as the eikonal limit of the standard Cutkosky rules. It is 
illustrated in figure 4 below. More explicitly, applying the prescription (2.16) to eq. (2.15), 
the imaginary part is determined as follows, 


2iImF (1) = (2ir) 2 ig 2 C F 9(v^)9{v%)(v 1 ■ v 2 )p 2e J ~ • ( 2 - 17 ) 

4 The momentum-space representation in eq. (2.15) is straightforwardly obtained from the eikonal Feyn¬ 
man rules. By Schwinger parametrizing the eikonal propagators in eq. (2.15) and performing the resulting 
Fourier transform, one recovers the position-space representation in eq. (2.1). 
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This representation of the imaginary part of the one-loop diagram motivates two remarks. 

The first remark concerns the region of momentum space which gives rise to the imag¬ 
inary part. Defining the light-cone variables k± = ^(fc° ± A: 3 ) and choosing the Lorentz 
frame in which the transverse components of the velocities vanish, ViT = 0, the support of 
the delta functions in eq. (2.17) is the region where the momentum of the exchanged gluon 
is maximally transverse, 

k T > k + ~ hr « 0 , ( 2 . 18 ) 

which was identified in ref. [21] as the Glauber region [70]. This agrees with the discussion 
in section 1: the imaginary part of eikonal diagrams arises from the exchanges of Glauber- 
region gluons. 

The second remark concerns the physical interpretation of applying the momentum- 
space cuts (2.16). By writing the delta functions in eq. (2.17) in terms of the plane-wave 
representation 5(A) = 2 J [f° d ue luA and performing the Fourier transform we find 

/ OO [* oo 

dti / dt 2 -(2.19) 

-oo J oo [— (tiVi - t,2V2 ) 2 + irj\ 

We observe that the resulting integration bounds compared to those of the uncut diagram 
in eq. (2.1) are extended according to dti —> dt*. This state of affairs can be simply 

understood on physical grounds: as the hard partons have been put on shell through the 
cutting rule (2.16), they are now asymptotic states propagating from fj = —oo to the 
interaction point. 



Figure 4. Graphical representation of cuts of eikonal diagrams in momentum space. 


We see that the position- and momentum-space representations of eikonal diagrams 
offer complementary points of view on the origin of their imaginary part. To summarize, in 
the position-space representation, the imaginary part is seen to arise from the exchanges 
of lightlike soft gauge bosons whose emission and absorption change the phases of the 
hard-parton states. In contrast, in momentum space, the imaginary part (related to the 
branch cut discontinuity through eq. (3.14)) arises from the two hard partons going on shell 
and exchanging Glauber gluons. Thus, the position- and momentum-space representations 
explain the origin of the imaginary part from the points of view of causality and unitarity, 
respectively. 

The momentum-space cutting prescription in eq. (2.16) has the conceptual advantage of 
factoring eikonal diagrams into on-shell lower-loop and tree diagrams which in turn can be 
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computed as independent objects. However, the resulting cut diagrams involve integrations 
over two-, three-, four-, ... particle phase space, as illustrated in figure 5. In practice, the 
evaluation of these phase-space integrals poses a substantial computational challenge which 
limits the applicability of the cut prescription (2.16) for obtaining imaginary parts. 


2 i Im 






Figure 5. The non-vanishing momentum-space cuts of a non-planar three-loop ladder diagram. 
The cuts require the evaluation of two-, three- and four-particle phase-space integrals. 


As we shall see in section 3, in position space, eikonal diagrams without internal ver¬ 
tices take the form of iterated integrals. In this representation, their imaginary parts can 
therefore be straightforwardly obtained by applying the principal-value formula ( 2 . 10 ) re¬ 
cursively. 

3 Position-space cuts of eikonal diagrams without internal vertices 

For completeness, in this section we review the derivation presented in ref. [1] of the imag¬ 
inary part of L-loop eikonal diagrams without internal (i.e., three- or four-gluon) vertices 
to the leading order in the dimensional regulator e. We will interchangeably refer to these 
diagrams as ladder-type diagrams. The basic observation is that in position space these 
diagrams are iterated integrals, and as a result their imaginary part can be obtained by 
decomposing the real-line integrations into principal-value and delta function contributions. 

In position space, an arbitrary L-loop eikonal diagram without internal vertices is 
composed of L soft-gluon propagators, interchangeably referred to here as rungs. Each 
rung extends between the Wilson lines spanned by any two (possibly identical) external 
four-velocities v\,...,v n where 1 < n < L + 1. For the jth rung we will denote these 
four-velocities by V(_. and v r . . We let denote the position of the kth attachment on the 
Wilson line spanned by Vi, counting from the hard interaction vertex and outwards, so that 
0 < ti t i < ti t 2 < < ti t Ni -1 < U,Ni> where N{ denotes the total number of soft-gluon 

attachments on the Wilson line. In addition, for the jth rung, we let the variables nij and 
rij record the soft-gluon attachment numbers on the Wilson lines spanned by V£. and v Tj , 
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respectively. The L-loop eikonal diagram is then defined as the 2L-fold iterated integral 


F {l) = C^ L \gp e ) 2L \\ 

3 =1 


d tg- mjdt r 


( v Cj ' v rj) rii=l rifc=0 ti,k) 

[—(.UjwVij - t rjtnj v rj ) 2 + irj] 


(3.1) 


where the kinematics-independent prefactor is determined by the color structure of the 
diagram and where it is implied that ij^+i = oo and ti t o = 0. Without loss of generality, 
we will assume that any rungs with both endpoints attached to the same Wilson line have 
been integrated out, and we suppress the resulting pole factors in e. (The additional factors 
produced by the integrations, involving epsilonic powers of the remaining variables, will not 
be of importance here, as our aim is to extract the imaginary part of F^ to the leading 
order in e.) 

To extract the imaginary part of F^ L > from the integral representation in eq. (3.1) it 
turns out to be useful to perform a change of variables which leaves each soft propagator 
dependent on a single variable. To this end, we adopt a change of variables introduced in 
ref. [11]. The idea is to first express the attachment points of the jth rung in terms of 
“polar” coordinates measuring the distance pj to the cusp (in units of the infrared cutoff 
1/A) and xj essentially measuring the emission angle of the soft gluon to the Wilson line, 5 



where 


0 < Pj < oo 
0 < x j < 1. 


After this change of variables, the diagram takes the form 


f (F = C iL \gp e ) 2L Yl 


ro ° d Pj f 1 

Jo p) 2t Jo 




3 = 3 -- '3 

where the soft propagators are defined as 


P%(?) 


Vi ■ Vj 


[—{xvi - (1 - x)vj) 2 + irj] 1 


and where the nesting of the integrations is encoded in 0, defined through 


(3.2) 


(3.3) 


(3.4) 


n Ni 

0(p, x) = 9(t i>k +i - t i}k ) 

i=l fc=0 



(3.5) 


We observe that the soft propagators’ dependence on the radial coordinates pj has scaled 
out in eq. (3.3), and that each propagator now depends only on a single variable Xj. This 
turns out to be particularly advantageous for the purpose of extracting the imaginary part 
of the diagram, as this circumvents the need to divide a higher-dimensional domain of 
integration into subdomains characterized by supporting a specific number of propagator 
roots. 

5 For a given rung, the two endpoints may of course be referred to interchangeably as left or right. 
However, for practical calculations, one particular choice may prove slightly more convenient. We refer to 
section 4 for examples. 
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Now we extract the overall infrared divergence of the diagram by setting n = p\ and 
then applying the following sequence of L — 1 substitutions 




tl 


Ul-i 

~ VL-l 


with 


0 < Tj < oo 
0 < Vj < 1 


where the variables tj have the dimension of length and the y.j are dimensionless. The 
L-loop eikonal diagram then becomes 



dxj p ef rj ( x i) K(x 1 ,...,x L ), 


(3.7) 


where the infrared divergence of the diagram is now absorbed into the kernel 


K(x i, ...,x L ) = g 


2 L 


n^) lu 


L—l 

n 

3 =1 ' 


dyjVj 


-1+2 je 


(1 - Vj) 1+2e @{{y,x}) . (3.8) 


Here 0({y, x }) denotes the result of applying the substitutions (3.6) to eq. (3.5). In analogy 
with section 2, we have here regulated the infrared divergence in a gauge invariant way 
through the exponential damping factor e~^ TL with A< 1. Eq. (3.8) contains in addition 
any potential ultraviolet subdivergences of the diagram (generated by the nesting function 
© ({!/.*}))• 

Having brought the L-loop eikonal diagram in the form (3.7), we now turn to extract¬ 
ing its imaginary part. Restricting our attention to the leading order in the dimensional 
regulator e, we will drop the dependence of the soft propagators on e, 



dxj P £% (xj)K(xi, ...,x L ) + 0(e d+1 ), 


(3.9) 


where d denotes the degree of divergence of the diagram, F ^ ~ T x (finite). 

To compute the imaginary part of eq. (3.9), we start by observing that eq. (3.8) is 
manifestly purely real. As a result, the Feynman ir/’s are the only source of imaginary 
parts of Eq. (3.9). Each of the x^-integration paths can therefore be decomposed into a 
principal-value part and small semicircles around the propagator poles. Given that the 
integrand takes purely imaginary values in the regions close to the poles and is real-valued 
on the remaining domain of integration, the resulting 2 L terms (which each involve L 
integrations) will be either purely real or purely imaginary. 

To collect the imaginary contributions, we define the cut propagator 


A ij(x) = -7 TVi- Vj6((xVi - (1 - x)vj) 2 ) , 


(3.10) 


and furthermore p-fold cutting operator 


n fl p „! 

G Ut. 2;, i .... i P n PV / dxj P(xj) / dx ik A(x ik ) K(xi,... ,xl) ■ (3.11) 

i =i k=1 , ' 0 


j+i.-.'p 
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The action of this operator is to replace the p propagators that depend on the specified 
variables by delta functions and to place a principal-value prescription on the integrals 
over the remaining variables. To simplify notation, we here dropped the indices on the 
(cut) propagators: P(xj ) = Pf\.{xj) and A (xj) = A .g r (xj). 

The imaginary part of any L-loop eikonal diagram without internal vertices can then 
be written, to the leading order in e, 

L L 

ImF (L) = E i p ~ l Cut,,,. F (l '> . (3.12) 

P= 1 ii,...,i p =l 

P odd 

This is the central formula of our approach [1], The formula (3.12) is illustrated schemati¬ 
cally for a generic ladder diagram in figure 6 below. 

We note that the decomposition of the line integrations in eq. (3.9) into principal- 
value and delta function contributions immediately shows that the imaginary part of the 
integrated expression for the eikonal diagram will have transcendentality weight one less 
than the real part. This follows from the fact that the delta functions will map the rational 
integrand to a rational expression after being integrated out. Thus, compared to the real- 
part contribution with L principal-value integrals, the weight is dropped by one. 



Figure 6. Schematic illustration of the formula (3.12) for the imaginary part of an eikonal diagram 
without internal vertices. The black dots at the endpoints of a soft-gluon propagator denote that the 
propagator has been cut—that is, replaced by a delta function. It is implied that the integrals over 
the attachment points of uncut soft propagators are principal-value integrals. The relative signs of 
the diagrams are determined by the factor i p_1 ; each individual diagram displayed here corresponds 
to the action of the cutting operator (3.11) on the eikonal diagram. 


It is natural to ask about the relation of the imaginary part of the eikonal diagram to 
the discontinuities in its various kinematic channels. This in turn leads us to ask for an 
appropriate set of variables in terms of which to express integrated results. A good choice 
of variables turns out to be given by the exponentials of the cusp angles, 

Xij = e~ lij , (3-13) 

with the cusp angles defined through cosh 7 ^ = | V{ ■ vj |. Expressed in terms of the Xij~ 
variables, the eikonal diagram has branch cuts located on the real line and satisfies Schwarz 
reflection, F^ L \xij) = P^ L Hxij)- As a result, the discontinuities of the diagram give rise 
to the imaginary part through the relation 

L 

2ilmF^ L \x) = E ' Vr i ) Disc i?(L) (*) • (3.14) 

1 \£ava 
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Here, the step functions account for the fact that the imaginary part has vanishing contri¬ 
butions from channels with space-like kinematics V( :j ■ v r . < 0. (This follows from the fact 
that propagators stretched between mutually space-like eikonal lines have vanishing cuts, 
as will be explained below eq. (3.16).) We will see an explicit example of this in section 4.3 
where we study a diagram that depends on two distinct cusp angles in purely time-like as 
well as mixed time- and space-like kinematics. 


In section 4 we will work out examples of how eq. (3.12) is used in practice to compute the 
imaginary part of ladder-type eikonal diagrams. To this end it will be useful to record the 
following partial-fractioned expressions, setting x = Xiji 




R(x) 


l 


l 


x — pi + ip x — f )2 — ip 


A 


(x) = ( S ( X ~ Pi) + S ( x - A*)) » 




(3.15) 


where the prefactor is the rational expression R(x) = x _ 2 = cothy^, and the denominator 
roots are given by 


(P1,P2) = < 


( -X _ O 

VX+1’ X+1J 


_x _L_ 

x~i ’1-x 


for Vi ■ Vj = cosh 7 ij > 0 

for Vi ■ vj = — cosh 7 ^ < 0 . 


(3.16) 


We note that in the upper case of eq. (3.16), the roots satisfy 0<p\<^<p2<l, whereas 
in the lower case they satisfy p\ < 0 < 1 < p 2 - Since the delta functions in eq. (3.11) are 
integrated over the interval [ 0 , 1 ], we may thus infer that the eikonal diagram will only have 
contributions to its imaginary part from channels with time-like kinematics Vi ■ Vj > 0 , as 
encoded in eq. (3.14). This is in agreement with the causality considerations of section 2. 

In section 4 we will make extensive use of the fact that the result for an eikonal diagram 
in time-like kinematics can be immediately obtained from the space-like result by analytic 
continuation of the cusp angle. To see this, let us first recall that the soft propagator takes 
the same form (3.4) in space- and time-like kinematics when expressed in terms of Vi ■ Vj , 
owing to our convention that all velocity vectors are outgoing. However, once expressed in 
terms of the relative angle 7 = 7 ^, it takes the respective forms 

— cosh 7 [— x 2 — (1 — x) 2 — 2 x(l — x) cosh 7 )] 1+e foru* 

cosh 7 [—x 2 — (1 — x ) 2 + 2 x(l — x) cosh 7 + ?r/] 1+e forxj 


P : 


- 


1-3 


X = 


Vj = — cosh 7 < 0 
Vj = cosh 7 > 0 , 


(3.17) 

where we dropped the ip in the space-like case, as the propagator roots are located outside 
the range [0,1] of x. Comparison of these expressions shows that we can map space-like to 
time-like kinematics by means of the analytic continuation 


— cosh 7 —> cosh 7 + ip , 
or equivalently, in terms of x = e -7 , 

1 

X —> - ip- 

X 


(3.18) 

(3.19) 


- 15 - 








4 Examples 

The aim of this section is to apply the formalism reviewed in section 3 to compute the 
imaginary part of a number of ladder-type eikonal diagrams. The main point to be ad¬ 
dressed here concerns the evaluation of the principal-value integrals involved in the p-fold 
cuts in eq. (3.11). 

4.1 The non-planar two-loop ladder diagram 

As a first example we will consider the non-planar two-loop ladder diagram, illustrated in 
figure 7 below. This diagram contains no ultraviolet subdivergence and therefore only has 
a simple pole in the dimensional regulator e. In agreement with the observations at the 
end of section 3, the diagram will only have an imaginary part for tinre-like kinematics 
v i • v 2 > 0. We therefore restrict our attention to this case. Since the diagram contains only 
one cusp angle, we will drop the subscripts for convenience and define coshy = v\ ■ v 2 as 
well as x = e -7 . 



Figure 7. The non-planar two-loop ladder diagram. 


The non-planar two-loop ladder diagram has the position-space representation 

d*i,i d* 1)2 df 2 .i d* 2)2 0 (*i ,2 - * 1 , 1 ) 0 (* 2,2 - t 2 ,i) (vi ■ v 2 ) 2 


F (2) = 


f 


[-(*i,Wi - t 2:2 v 2 ) 2 + irj] 1 ' [-(fi, 2 t>i - t 2) iv 2 ) 2 + ir /] 


1—e ’ 


(4.1) 


where the prefactor is given by C^ = — 3 2 < n F P ■ To compute the imaginary part of 

this diagram, our first task is to write it in the form of eq. (3.7). This is achieved through 
the change of variables in eq. (3.2), followed by that in eq. (3.6), 


ii,i 

*2,2 


= A 


X 

1 — X 


* 1,2 

*2,1 


= a 


y 

i -y 


followed by 


A 


= P 


S- t 


After these transformations the diagram takes the desired form, 

F (2) = f dx I' d y Py}(x) P®(y) K(x, y). 

Jo Jo 


(4.2) 


(4.3) 
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where the kernel K(x,y), upon the additional change of variable u = is given by 


K(x,y) = y, 4e J^ d ^ e _ 4e ^ duu 1+2e (u + l) 4 e d(^-uj 0 (u - (j —. (4.4) 

By comparing eqs. (4.1) and (4.3), we see that the effect of the first two transformations 
of eq. (4.2) is to leave each soft propagator dependent on a single variable. The effect of 
the last change of variable is to extract the overall infrared divergence of the diagram. 

To facilitate the computation of the cuts in eq. (3.12) we will first evaluate the integral 
K(x,y). The u-integral in eq. (4.4) may be performed in terms of the 2 -Pi hypergeometric 
function. The primitive has the e-expansion 

f(u) = - 7^7 2 -Pi( 2 e, 4e; 1 + 2e; —u) = (l + 0(e 2 )) = ^ + log u + 0{e), (4.5) 

and so K(x,y ) has the e-expansion 


K(x,y) = T(4e) (£) 6(y - x) 


/(-)-/^ 
XJ \ 1 — X 


1 

4e 


(x) ~ X ^ ( log x ~ log ir|) + ' ( 4 - 6 ) 


Substituting this result for K(x,y) into eq. (4.3), we can write the non-planar two-loop 
ladder diagram in the convenient form 


p ( 2 ) = 



(4.7) 


where F^ is finite, given to leading order in e as, 

n 1 pi 

= J d x J d y (x) P$ (y) (d(y - x) - 0{x - y)j log |. (4.8) 

Here we have dropped the dependence of the soft propagators on e and furthermore rewrit¬ 
ten the integrand of F O to make its symmetry under x <—> y manifest. (This was achieved 
by changing variables (x,y) 1 —> (1 — x, 1 — y) on the second term log arising from 
eq. (4.6).) 

As the prefactor of F ^ in eq. (4.7) is real, it factors out on both sides of eq. (3.12), 
yielding the formula 


Im F® = Cut x F {2) + Cut y P (2) , 


(4.9) 


where we recall that Cut Xi is defined in eq. (3.11) and replaces the propagator depending 
on the specified variable by a delta function and places a principal-value prescription on 
the integral over the remaining variable. The two cuts on the right-hand side are equal 
because the integrand of eq. (4.8) is symmetric under the interchange of x and y. (This 
also follows from the v\ <—> V 2 symmetry of the original diagram.) Thus, it suffices to 
compute Cut x F^ 2 \ given by 

pi pi 

Cut X F (2) = j dxPV J dy A 12 (x)P$(y) (d(y - x) - 6(x - yj^log^ . (4.10) 
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Inserting the partial-fractioned expressions for Ai 2 (x) and P[ 2 (y) given in eq. (3.15) and 
performing the trivial integral over x produces 


Cut x FW = ~R(x) 


' 2 E 

k =1 


PV 


d y - 


rpk 


d y 


1 pk 


i 


y-pi+ip y-p 2 -ir] 


)log —, (4.11) 

> Pk 


where the propagator roots p \.2 are given in the upper part of eq. (3.16). 

We are now confronted with the task of evaluating principal-value integrals. As such 
integrals do not immediately take the form of iterated integrals, our strategy for evaluation 
will be to write them as differences of iterated integrals, which in turn are readily express¬ 
ible in terms of multiple poly logarithms. The basic observation is that the principal-value 
integral equals the corresponding full integral minus the imaginary part of the latter, cf. 
eq. ( 2 . 10 ). 

As a simple illustration, let us consider the evaluation of the following principal-value 
integral, 


PV 



d y 

y- pi+irj 


f 1 d y 

o y-pi + iv 



d y 

y - pi + ir) ' 


(4.12) 


The full integral evaluates to G(pi\ 1) by definition, cf. eqs. (A.1)-(A.2). Its imaginary part 
arises from the pole of the integrand and is extracted by localizing the integration variable, 


Im 


d y 


I o y-pi+ir) 


= —7r 


d yS(y- pi) = - 7 r. 


(4.13) 


where we used in the last step that the pole is located inside the range of integration, in 
agreement with the discussion below eq. (3.16). Thus we arrive at 


PV T- ^— = G(pi-l) + 7Ti. (4.14) 

Jo y- pi+ 

In this simple example we could have computed the real part more directly, 

pv ^ ^ = R eG(pr; D = Itelog (l - £) = log |l - £| = lo g (£ - l) . (4.15) 

However, an extension of this direct approach to higher-weight cases requires the use of a 
sequence of functional identities which in practice is case-dependent and thus not applicable 
in a systematic way. In contrast, the above method relies only on a construction of the 
imaginary part which can be derived systematically as demonstrated in appendix A. 

Returning to eq. (4.11) and evaluating the principal-value integrals following the steps 
outlined above, the result for the cut is readily expressed in terms of multiple polyloga¬ 
rithms, 

Cut* PW = ~R( X ) 2 (~2G(p2, 0; 1) + 2G( P1 , 0 ; 1 ) - 2G(0, p 2 ; P 2 ) - 2G(0, p 2 \Pi) 

+ 2G(0, pi; p 2 ) + 2G(0, pi\pi) + G(0; p 2 )G(p2\ 1) + G(0; pi)G(p 2 ] 1) 
- G(0; P2 )G{pi- 1 ) - G(0; pi)G(pi; 1) - 2vriG(0; |)) , (4.16) 
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where we refer to eqs. (A.1)-(A.2) for definitions. In this expression, the multiple poly¬ 
logarithms depend on the propagator roots Pk(x) through both their indices and their 
arguments. This expression can in turn be rewritten in terms of polylogarithms with con¬ 
stant indices by exploiting the Hopf algebra structure of multiple polylogarithms, which 
encodes the plethora of functional identities within this class of functions [71-76]. Utilizing 
this algebraic structure, we have implemented the steps required to achieve the desired 
functional form as a general algorithm. We refer to appendix B for further details. The al¬ 
gorithm leaves a simplified form of eq. (4.16) expressed in terms of harmonic polylogarithms 
which, using eq. (A. 5), can be simplified further into classical polylogarithms, 

Cut, F™ = -vr R( X ) 2 (~G(0 ,1; X ) + G(0,0 ; x ) - G(0, -1; x) - ^ 2 ) 

= -7^(x) 2 (Li2(x 2 ) + l°g 2 X- C 2 ) ■ (4.17) 

With this result, we can now immediately obtain the imaginary part of the two-loop ladder 
from eq. (4.9), recalling that the two cuts are equal. We find 

ImT (2) = -7ri?(x) 2 (Li 2 (x 2 ) + log 2 x - C 2 ) , (4.18) 

and we recall that multiplying the infrared pole cf. eq. (4.7) onto both sides of eq. (4.18) 
gives the imaginary part of the original diagram F ( 2 ). This completes the evaluation of the 
imaginary part of the non-planar two-loop ladder diagram to the leading order in e. 

As a crosscheck of the result in eq. (4.18), we can alternatively compute the imaginary 
part of the non-planar two-loop ladder by evaluating the diagram for space-like kinematics 
v\ • V 2 < 0, in which case it will be purely real (cf. the discussion at the end of section 3), 
and subsequently perform the analytic continuation to time-like kinematics. We refer to 
the end of section 3 for a more detailed discussion of analytic continuations. 

To the leading order in e, the two-loop ladder is given by eq. (4.8), although we must 
bear in mind that for space-like kinematics the propagator roots pk are given by the lower 
case of eq. (3.16). Inserting into eq. (4.8) the expressions for pf!} given in eq. (3.15), the 
diagram readily evaluates into multiple polylogarithms, 

^(2) = ( G(/0l; 0 , Pl -1 ) - G( pil 0, p 2 ; 1) - G(p 2 , 0, Pi; 1) + G(p 2 , 0, p 2 ; 1)) . (4.19) 

The tilde on the left-hand side serves to remind us that the expression for the diagram 
on the right-hand side is valid for space-like kinematics. We can use the algorithm in 
appendix B to recast this representation in terms of polylogarithms with constant indices. 
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In fact, the two-loop ladder diagram can be expressed in terms of classical polylogarithms, 6 

R {2) = R(X? ( 2 G( 0 , 1 , 0 ; x) ~ 2G(0, 0 , 0 ; x ) + 2G(0, -1,0; x) - C 2 G'( 0 ; x ) - Cs) 

= R(X) 2 (Li.3(x 2 ) - log X Li 2 (x 2 ) - g log 3 x - C '2 log X - ( 3 ) ■ (4.20) 

We can now find the result for the two-loop ladder diagram in time-like kinematics by 
performing the analytic continuation x — > — 1/x — ' l/r l ° n eq. (4.20). Under the analytic 
continuation, the rational function R(x) picks up a minus sign, while the polylogarithms 
transform according to 

log X —> - log X - ni 

Li 2 (x 2 ) - Li 2 (x 2 ) - 2 log 2 x + 2C 2 - 2m logx 

Li 3 (x 2 ) Li 3 (x 2 ) + | log 3 x - 4C 2 log X + 2m log 2 x • (4.21) 

Applying these replacements to eq. (4.20) we find the following result for the non-planar 
two-loop ladder with time-like kinematics, 

R {2) = R(x? ( Li 3 (x 2 ) - log X Li 2 (x 2 ) - l log 3 X + 5C 2 logX - C 3 

- ivr(Li 2 (x 2 ) + log 2 x - C 2 )) • (4-22) 

We observe that the imaginary part of eq. (4.22) agrees with the result found in eq. (4.18), 
as expected. We conclude that the cutting prescription for the two-loop ladder stated in 
eq. (4.9) produces the correct imaginary part. The cutting prescription (4.9) is illustrated 
in figure 8. 



Figure 8. Graphical representation of the cutting prescription for the non-planar two-loop ladder 
stated in eq. (4.9). The black dots at the endpoints of a soft-gluon propagator indicate that the 
propagator has been cut; i.e., replaced by a delta function. It is implied that the integrals over the 
attachment points of uncut soft propagators are principal-value integrals. 


6 Note that in the representation in the second line of eq. (4.20), the polylogarithms and logarithms have 
the respective arguments y 2 and This form has the advantage of being particularly compact as well as 
convenient for the purpose of performing the analytic continuation \ — 1/x — which then maps the 

arguments of the (poly)logarithms to the vicinity of the branch cuts. (In general, computer algebra software, 
such as Mathematica augmented with the package HPL [77, 78], does not detect the 2 7r monodromy of log(y 2 ) 
after tracing out a complete circle around the branch point.) 
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4.2 Three-loop non-planar ladder diagram 

To demonstrate that the principal-value integrals involved in the p-fold cuts in eq. (3.11) 
can indeed be evaluated in non-trivial cases, we consider in this section the three-loop 
ladder diagram illustrated in figure 9. This diagram also represents an example of an 
eikonal diagram with multiple-cut contributions to its imaginary part (in the case at hand, 
a triple cut). As in section 4.1, we take coshy = v\ -V 2 > 0, in order to have a non-vanishing 
imaginary part, and set x = e_7 - 



Figure 9. A three-loop non-planar ladder diagram. 


The position-space representation of the diagram in figure 9 takes the form 


p(3) 



nn dt 

. i= 1 j =1 


UiJ=l +1 - tjj) (ui • V 2 ) 3 

D{t\ t 1 , t 2 , 3 ) ^2,l) ’ 


(4.23) 


where D(t\,t 2 ) = [— (t\vi — £ 2 ^ 2) 2 + irj ] 1 c . To compute the imaginary part of this diagram 
from eq. (3.12), our first task is to bring it into the form of eq. (3.7). This is achieved 
through the changes of variables in eq. (3.2) (with t£ jjTnj = tij), followed by the sequence of 
substitutions in eq. (3.6), setting (aq, aq, £ 3 ) = (x, y, z) and (yi, 2 / 2 ) = (t, u ) for convenience. 

After these transformations, the diagram takes the form 


f (3 ) = c (3) f dx dy dz p H ^ pH ^ ( z )K(x,y,z ), 


(4.24) 


where the kernel is given by 


K(x,y,z) = T(6e) J dt t~ 1+2e (l - t)~ 1+2e j da T 1+4e (1 - it) 


—l+2e 


1 y 1—u 


a I y _ t_ \ nl tu _ l-z \ nl 1-z _ (l-t)u 

u lx l-x) u \l-y l-i 


(4.25) 


The arguments of the step functions simplify after rescaling the integration variables ac¬ 
cording to p-j: 1 —y t and 1 —> u. As we are interested in computing the imaginary part of 
F^ only to the leading order in e, we may set e to zero in the u-integral. Performing the 
t-integral in terms of hypergeometric functions and subsequently expanding in e, we find 
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the expression 

A-( w ) = i(|) 

_ 1 /^\ 6e 


6{y — z)6(z — x)\ log 2 
JC(x,y,z), 


1 — X z 
x 1 — z 


+ ( y <—> z) 


6 eVA J - ( 4 ‘ 26 ) 

valid to the leading order in e. We observe that eq. (4.26) is symmetric under the interchange 
of y and z; hence the integrand of the full diagram in eq. (4.24) is as well. (In other 
words, interchanging the two parallel gluon lines in figure 9 leaves the diagram invariant.) 
This observation implies that Cut^ F ® = Cut 2 and thereby reduces the number of 
independent cuts to be computed. 

Substituting eq. (4.26) into eq. (4.24) we can write the non-planar three-loop ladder 
diagram in the convenient form 


_ 


where pC) is finite, given to leading order in e as 


F (3) = dxdy dz (®) P S (v) P S 0 ) V, z ) ■ 


(4.28) 


As the prefactor of pC) in eq. (4.27) is real, it factors out on both sides of eq. (3.12), 
yielding the formula 

Im F {3) = Cut* P (3) + Cut y F {3) + Cut z P (3) - Cut x ^ z P (3) . (4.29) 


The appearance of a multiple-cut contribution represents a new feature for diagrams with 
more than two loops. Incidentally, the diagram F^ is the only three-loop ladder diagram 
with an O(j) divergence that has a non-vanishing triple cut. This point is illustrated in 
figure 10. The diagram F^ therefore provides an excellent example to demonstrate a 
multiple-cut contribution. 

Let us start by evaluating this triple cut: as all three integrations are localized by delta 
functions, the cut is immediately computed, 

Cut x ,y,zF^ = f dx dy dz Ai 2 (x) A 12 (y) A i 2 (z) )C(x, y, z) 

Jo 

3 2 

= ~Y P (x) 3 ^2 lC{pk,pi,p m ) ■ (4.30) 

k,l,m =1 

In the second line we inserted the form of A 12 given in eq. (3.15) and then integrated out the 
delta functions. Now observe that lC(pk, Ph Pm), which is implicitly defined in eq. (4.26), 
is non-zero if and only if ( k,l,m ) = (1,2,2). Indeed, in the first term of eq. (4.26), the 
logarithm is non-zero only when k 7 ^ m, while the step functions dictate that p/~ < p m < p\. 
Since pi < P 2 , we must therefore have k = 1 and l = m = 2 . An identical argument applies 
to the second term of eq. (4.26). We arrive at the simple result 

Cut*,^ P (3) = - yP(x ) 3 log 2 x • (4.31) 
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Vi V 2 



Vi V 2 



(b) 


Figure 10. Triple cut of the two three-loop diagrams that have a single pole divergence. In (a) the 
non-vanishing triple cut of F^ is illustrated by the exchange of three lightlike gluons. In contrast, 
(b) shows graphically that the triple cut of the maximally crossed diagram vanishes, because the 
three gluons cannot simultaneously be aligned with the light cone (i.e., go on shell). 


This completes the evaluation of the triple-cut contribution to eq. (4.29). 

We now turn to the single cuts. We will focus on the contribution Cut^pl 3 ', the 
remaining cuts being computed completely analogously. (Also recall from the discussion 
below eq. (4.26) that Cut^ J 7 ^ 3 ) = Cut z so that only one additional single cut needs to 
be computed.) Using the partial-fractioned expressions for A 12 and i n eq. (3.15) and 
integrating out the delta function we find, dropping ±ir; for notational convenience, 


Cut x J 7 ^ = /' dx PV f 1 dy PV I' dz A 12 (x)PS(y)PS(z)JC(x,y,z) 
Jo Jo Jo 

^<x) 3 E pv f M ^ 


7 r 


k= 1 


y-p 1 y~P2 


PV / dz 

J pk 


1 _1 

2-Pl 2 —P2 


(4.32) 


exploiting the residual y <—> z symmetry of the integrand. 

The principal-value integrals may be evaluated as the corresponding full integrals minus 
the imaginary part of the latter. Formulas of imaginary parts of multiple polylogarithms are 
listed up to weight four in appendix A. The first PV integral gives rise to step functions 
involving the variables pi,p 2 and y. Splitting the y-integral in the k = 1 term into two 
integrals with the respective domains [pi,p 2 ] and [/ 02 > 1 ] allows all step functions to be 
resolved. Ultimately, the cut is found to evaluate into the expression 


Cut* P (3) 



g l°g 4 X + log 2 x(h 2 (x 2 ) - 3 C 2 ) - 2 log x[h 3 (x 2 ) - C 3 ) 


+ h 2 ,2(x 2 ) + 2 h 4 ( x 2 ) ~ ( 2 h 2 ( x 2 ) ~ iC4 


(4.33) 
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This result is expressed in terms of harmonic polylogarithms, defined through eq. (A. 6 ). 
Similarly, the y- and z-cuts are found to take the form 

Cuty = Cut* J r(3) = jR{x) 3 -5 log 4 X + log 2 x(Ho(x 2 ) + K 2 ) 

+ 2F 3 ,i(x 2 ) + ^2, 2 (X 2 ) - C2^ 2 (X 2 ) + fC 4 I • (4-34) 


Combining all single and triple cuts according to eq. (4.29) yields the following imaginary 
part 


Im J-4 3 ) 


nR(xf 


-|l°g 4 x + l°g 2 x(^ 2 (x 2 ) + 3 C 2 ) - logx(^3(x 2 ) - C.3) 
+ i?3,i(x 2 ) + tf 2)2 (x 2 ) + ^(x 2 ) - C2^ 2 (x 2 ) + 2 C 4 • 


(4.35) 


As a crosscheck of this result, we can alternatively compute the imaginary part of the 
three-loop ladder by evaluating the diagram for space-like kinematics v\ ■ v 2 < 0 , in which 
case it will be purely real (cf. the discussion at the end of section 3), and subsequently 
perform the analytic continuation to time-like kinematics. 

To the leading order in e, the three-loop ladder is given by eq. (4.28), although we must 
bear in mind that for space-like kinematics the propagator roots pk are given by the lower 
case of eq. (3.16). Inserting into eq. (4.28) the expressions for p) ( | and }C(x,y,z) given in 
eqs. (3.15) and (4.26), respectively, the diagram is directly expressible in terms of multiple 
poly logarithms, 

^(3) = R^f £ (-iy+i+ k + l +™G(p i+1 ,p j+1 ,k,l,p m+1 -l), (4.36) 

i,j,k,l,m= 0,1 


where the tilde on the left-hand side indicates that the diagram is computed for space¬ 
like kinematics. We can use the algorithm in appendix B to recast this representation in 
terms of poly logarithms with constant indices. In fact, the three-loop ladder diagram can 
be expressed in terms of harmonic poly logarithms, 

R(x ) 3 


jr(3) = 


T 5 log 5 X + | log 3 x(# 2 (x 2 ) - C2) - 2 log 2 x(#3(x 2 ) - C 3 ) 

+ 4 log x (#3,1 (X 2 ) + tf 2 , 2 (x 2 ) + #4(X 2 ) + C 2 # 2 (X 2 ) + |C 4 ) 

- 6P 4 ,i(x 2 ) - 6 . 03,2 (x 2 ) - 4-0 2 ,3(x 2 ) - 6-05(x 2 ) - 2C 2 H 3 ( X 2 ) 

+ 4(3 -0 2 (x 2 ) + 3(5 + 2 C 2 C 3 I • (4.37) 


We have cross-checked this expression with previous results in the literature, finding agree¬ 
ment. 1 

7 More specifically, adding the diagram in figure 9 to the maximally-crossed three-loop ladder, which we 
have computed by the same methods, we find agreement with eq. (A.l) of ref. [28] (in its published version; 
or alternatively eq. (73) of the corresponding arXiv e-print (v2)). Furthermore, the two color structures of 
the (3,3) web are linear combinations of these two diagrams, according to eq. (4.26) in ref. [39]. Inserting 
our results for the two diagrams, we recover the color structures in their eqs. (4.29) and (4.33) (where the 
basis functions are given explicitly in appendix A), thereby cross-checking our results for both diagrams 
individually, and in particular our result in eq. (4.37) for the three-loop non-planar ladder. 
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We can now find the result for the three-loop ladder diagram in time-like kinematics 
by performing the analytic continuation x — 1/x — ir) on eq. (4.37). Under the ana¬ 
lytic continuation, the rational function R( X ) picks up a minus sign, while polylogarithms 
transform according to 


log X ->• log(— 1/x - iij) = - log x - m , 

Ha(x 2 ) H s ( 1/x 2 + iri) • (4.38) 


Thus, all harmonic polylogarithms are evaluated slightly above the branch cut [1, oo). They 
were subsequently expressed in terms of H$(x 2 ) and l°g X using the Mathematica package 
HPL [77, 78]. In this way, we find the following result for the three-loop ladder with time-like 
kinematics, 


jt(3) 


R(x ) 3 

4 


TE lo S 5 X - | log 3 x[h 2 (x 2 ) + HC2) + 2 log 2 x(h 3 ( X 2 ) ~ C3) 

- 41ogx(i7 3 ,i(x 2 ) + H 2 , 2 (x 2 ) + H A { X 2 ) - 5 ( 2 H 2 ( X 2 ) - f C 4 ) 
+ 677 4 ,i(x 2 ) + 677 3 , 2 (X 2 ) + 4i7 2 , 3 (x 2 ) + 6 H 5 ( X 2 ) - 10( 2 H 3 ( X 2 ) 


- 3 Cs + IOC2C3 

+ 4vrlog 4 x - log 2 x(# 2 (x 2 ) + 3C 2 ) + logx(-ff 3 (x 2 ) - C 3 ) 

- TMx 2 ) - h 2 , 2 (x 2 ) - ^4(X 2 ) + C2Ho(x 2 ) - ia) • 


4Cs^2(x 2 ) 


(4.39) 


We observe that the imaginary part of eq. (4.39) agrees with the result found in eq. (4.35), 
as expected. We conclude that the cutting prescription for the three-loop ladder stated in 
eq. (4.29) produces the correct imaginary part. The cutting prescription (4.29) is illustrated 
in figure 6. 


4.3 Two-loop web with three Wilson lines 

The formalism of section 3 allows us to compute the imaginary part of eikonal diagrams 
only to the leading order in e. This appears to limit the applicability of the approach, but 
in practice a large class of diagrams have only simple poles in e, and the coefficient of the 
4 pole of the correlator of Wilson lines defines physical observables of interest, such as for 
example the cusp anomalous dimension (i.e., the anomalous dimension of the correlator 
of two Wilson lines, cf. eq. (2.8)). As it turns out, the cusp anomalous dimension can 
be expressed entirely in terms of diagrams with simple poles in e (once the diagrams are 
expressed in terms of the renormalized coupling gp). We observed this already for Abelian 
gauge theories in eq. (2.7), but the statement extends to the case of non-Abelian gauge 
theories as well. This owes to the non-Abelian exponentiation theorem for two Wilson 
lines [4, 5] which states that the two-line correlator can be written as the exponential of 
the sum of webs, defined as the subclass of diagrams which are eikonal-line two-particle 
irreducible. The webs appear in the exponent with their color prefactors appropriately 
modified to account for the color factors of the complete set of diagrams arising from 
expanding the exponential. 
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Over the past four years, non-Abelian exponentiation has been shown to generalize to 
correlators of an arbitrary number of Wilson lines [8-14], In this section we will study the 
interplay of this remarkable theorem with the formalism of section 3. 

The statement of non-Abelian exponentiation in the multi-line case requires a new 
classification of the set of diagrams which appear in the exponent. In the multi-line case, 
a web is defined as a collection of diagrams which are mutually related by permutations 
of the order of gluon attachments, acting on each Wilson line separately. As an example, 
consider the ( 1 , 2 , 1 ) web in figure 11 , whose imaginary part we will turn to shortly. 



Figure 11. The two diagrams constituting the (1,2,1) web. The antisymmetric combination of 
the two diagrams is not captured by exponentiation of the one-loop subdiagrams, and this linear 
combination of the diagrams, defining the web, appears in the exponent. Each separate diagram 
diverges as 0{\), but the leading divergences conspire to cancel in the web, leaving an 0(4) 
divergence. 


The labels (1,2,1) refer to the number of gluon attachments on each eikonal line, 
and the two diagrams are related by permutations of the attachment points, where the 
permutations act on each line separately. The individual diagrams in figure 11 are clearly 
not eikonal-line 2PI diagrams. However, the contributions of these diagrams are not entirely 
reproduced by exponentiation of their one-loop subdiagrams. Indeed, the sum of the one- 
loop subdiagrams spanned by the pairs of lines {ui,^} and {^ 2 ,^ 3 } will appear in the 
exponent, but upon expanding the exponential to second order, these subdiagrams only 
reproduce the symmetric linear combination of the two-loop diagrams. To compensate, the 
antisymmetric linear combination of the two-loop diagrams, which defines the web, must 
be added to the exponent as a contribution. Denoting the color and kinematical factors of 
the diagrams in figure 11 as C^b and J-^b, the web is accordingly defined as 

^£ 2 , 1 ) = (C a , C'b) 7 ) - C 'b)(-^a - Tb) • (4.40) 

We can now state the general non-Abelian exponentiation theorem [8-14]. Recalling the 
definition in eq. (1.1) of the Wilson line spanned by the four-velocity v^, the statement 
is that the correlator of an arbitrary number of Wilson lines is given as the exponential of 
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the sum of webs, 

= expf £ Cj'RifX (4.41) 

^ i£ {webs} ' 

where each web i contributes to the exponent through the color C) and kinematical T) 
factors of its constituent diagrams, weighted by means of the web mixing matrix R{ in 
analogy with eq. (4.40). The web mixing matrices can be computed systematically by 
means of the replica trick of statistical mechanics [ 8 , 79]. Among several properties they 
satisfy the zero-sum-row condition Rab = 0 which ensures that the symmetric linear 
combination of the constituent diagrams is projected out [ 10 ]. 

The mixing matrices satisfy an additional weighted zero-sum column condition which 
ensures that the leading divergence of the constituent diagrams of a web conspire to cancel 
when the diagrams are added. This is a general feature of webs, ultimately following from 
their renormalization properties [ 8 , 9, 11], leaving in many cases a web with an O(^) 
divergence. As a result, webs are particularly amenable to the formalism of section 3. It 
is here important to keep in mind that the cutting prescription should be applied to an 
entire web rather than its constituent diagrams separately, as the separate imaginary parts, 
computed to leading order in e, will cancel. 

To illustrate the procedure in detail, we turn to the web in figure 11 and compute its 
imaginary part. In analogy with sections 4.1 and 4.2, we take cosh 712 = v\ ■ V 2 > 0 and 
cosh 723 = V 2 ■ V 3 > 0 , in order to have a non-vanishing contribution to the imaginary part 
from both kinematical channels, and set x = e -712 and ip = e -723 . Our first task will be 
to show that the leading 0 {\) divergences of the individual diagrams conspire to cancel, 
leaving an O(^) divergence. We will then apply the cutting prescription of eqs. (3.10)-(3.12) 
directly to the web written in a form with a manifest 0 (~) divergence. 

Let us start by considering diagram 6 a in figure 11. Its kinematical factor is given by 


jr _ (j{2) r dfyi dt 2 ,i dt 2 ,2 df 3 ,i 0{t 2 , 2 - t 2 ,i) (t>i • v 2 ) {V2 • V 3 ) 
Jo [-(ii,m - t2,2V2 ) 2 + ir ]} 1 e [-(t2,iV2 - t 3 ,iv 3 ) 2 + irj ] 1 


(4.42) 


where C ^ contains coupling constants etc., but no color factor. In analogy with sec¬ 
tions 4.1 and 4.2 our first task is to write this expression in the form of eq. (3.7). This 
is achieved through the changes of variables in eq. (3.2) (with = £ 2 , 3 -}), setting 

(aq, X 2 ) = (x,y) for convenience. 

After these transformations, the kinematical factor takes the form 


J- a = f'dxdy Pg (x) Pg ( y ) K(x, y ), 

Jo 

where the kernel is given by 


,,, \ 4 e f°° dpidp 2 Qt N 

K(x,y) = y, j ^ 1 —-^- e 6{p l x-p 2 y) 


r(4e) 

2 e 



[u 2 e 2 i ? 1 (2e,4e; 1 + 2e; -u)] 


U =OO 

u=y/x 


(4.43) 


(4.44) 
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The result for K (x, y) was obtained by applying the substitution in eq. (3.6) and performing 
the remaining integrations in complete analogy with section 4.1. The second diagram of 
the web, J 7 , differs only in the step function which reads 0(p2V — Pix), changing the lower 
integration bound in eq. (4.44) from y/x to x/y. After expanding the gamma function and 
the hypergeometric function in e, we thus find the kinematical factor of the ( 1 , 2 , 1 ) web to 
take the form 


P 3 , Pb — 


C( 2 ) 

8^ 


(x) / dxd y p S( x ) p S(y) 


2e 


2e 


(1 + 0(e 2 )). (4.45) 


Upon expanding [ • • • ] in e, we observe that the leading poles of the separate diagrams 
cancel, leaving an 0(~) divergence, in agreement with the discussion above. 

Factoring out the remaining pole, we can write the web in the convenient form 


w (2) c a -c h cw ( p\^ m 


(4.46) 


( 2 ) 

where P { ] 2 ^ is finite and given to the leading order in e by 


r-l 


x 


^(?2 i) = / dx d y Pi 2 J (x) P^ (y) log - . 

v ’ ’ ’ Jo y 


(4.47) 


(o) 

As the prefactor of P (l 2 ^ in eq. (4.46) is real, it factors out on both sides of eq. (3.12), 


yielding the formula 


Im= Cut ^- 7 r (l 2 2 ,l) + Cut y '^""( 1 , 2 , 1 ) 


(4.48) 


More explicitly, by inserting the definition of the operator Cut Xi in eq. (3.11), we have 


III! P, 


(12 1 )= [ dxPV [ dy A u (x)P^{y)log- + PV [ dx [ dy P$(x) A 23 (y) log 
v ’ ’ ’ Jo Jo y Jo Jo 


y 

(4.49) 


We observe that the second term equals minus the first term with the two cusp angles and 
integration variables interchanged, making the imaginary part of the web antisymmetric 
under the interchange 712 <—> 723 - (The antisymmetry is of course inherited from the 
full web which has this property by construction.) This observation allows us to write the 
imaginary part in the manifestly antisymmetric form 

Im ^( 1 , 2 , 1 ) = J (vV0 - J(i>,x), (4.50) 

where the auxiliary function is defined as the first term of eq. (4.49), 

p 1 p\ 

J{x,i>)= l dx PV / dy A u (x)P^{y) log-. (4.51) 

Jo Jo y 
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Evaluation of the imaginary part of the web thus reduces to the evaluation of the inte¬ 
gral J(x, ip). The latter can be computed by recalling the partial-fractioned expressions for 
A 12 and P ijg given in eq. (3.15) which in the present notation read 

Ai 2 (x) = - pi(x)) + S(x - p 2 (x))) , 

^s 1 (2/) = ™ -wry-—) , (4-52) 

2 \y-pi(V’) + M 7 y - P 2 {v) - Hi J 

where we wrote out the expressions explicitly to emphasize their dependence on the two 
distinct kinematical invariants x and ip- 

Integrating out the delta functions in eq. (4.51) leaves one principal-value integral to be 
evaluated. This integral is computed as the corresponding full integral minus its imaginary 
part, as explained in section 4.1. 8 In this way we find 

J(X, VO = -jR( X )RW £ PV [ 1 d y ( -\- f . ) log ( Pkix) ) 

4 Jo \y-piw) + ir i y ~ «(w - ^) \ v ) 

= ~R(x)RW (-4Li 2 (-ip) + log 2 ip - 41og ip( logi'tp + 1) - log(x + 1)) 

- 2 log X log Ip ~ 2 C 2 ) • (4.53) 

Upon the antisymmetrization in eq. (4.50) the terms on the last line of eq. (4.53) cancel, 
and we find the following result for the imaginary part of the (1, 2,1) web, 

Im = _7r R (x)RW (Li 2 (—x) - Li 2 (—VO - \ (log 2 x - log 2 ip) 

+ (log X + log VO (l°g(x + 1) - log(V> + 1))) • (4.54) 

As a crosscheck of this result, we can alternatively compute the imaginary part of the web 
by evaluating the diagram for space-like kinematics v 1 ■ u 2 <0 and u 2 • V 3 < 0, in which case 
it will be purely real (cf. the discussion at the end of section 3), and subsequently perform 
the analytic continuation to time-like kinematics. 

To the leading order in e, the web is given by eq. (4.47), although we must bear in 
mind that for space-like kinematics the propagator roots pk are given by the lower case of 
eq. (3.16). We now insert into eq. (4.47) the expressions for given in eq. (3.15) and 
perform the integrals in analogy with the calculations in eqs. (4.19)—(4.20) for the case of 
the non-planar two-loop ladder. This leads to the following result for the web in space-like 
kinematics, 

^( 1 2) 2 , 1 ) = R(x)RW {Hip) log x - Hx) log VO , (4.55) 

where we introduced the auxiliary function 

Hx) = ~ Li 2 (l - x) “ il°g 2 X- (4-56) 

8 To proceed we assume, without loss of generality, that ijj < This fixes pi(x) < Pi (VO < P2(V0 < P2(x)> 
which allows the step functions in eqs. (A.22)-(A.25) to be resolved. 
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These expressions are consistent with results previously obtained in the literature, see for 
example eq. (3.11) in ref. [38], as well as references therein. We can now obtain the result 
for the web in time-like kinematics by performing the analytic continuations x — 1 /x ~ iv 
and i)j —> — l/ip — irj on eq. (4.55). Under the analytic continuation, the functions appearing 
in eq. (4.55) transform as 


R(z) -> -R{z) 
log z —> — log z — ni 

L(z) -► -Li 2 (-z) - log z log(z + 1) + | log 2 z - 5 C 2 

- 7 ri( log(z + 1) - 2 logz) • (4-57) 


Upon analytic continuation in x and V’ we thus find the following result for the web with 
time-like kinematics, 


77(2) 

‘'( 1 , 2 , 1 ) 


R( X )RW 


- wr(Li 2 (—x) - \ log 2 X + (log X + log V’) log(x + 1 )) 
+ log ,\Li 2 (—VO - (log X log ip - 6C2) log(x + 1 ) 

+ | (log X log ip - IOC2) log x - (x <—^ 4 >) ■ 


(4.58) 


We observe that the imaginary part of eq. (4.58) agrees with the result found in eq. (4.54), 
as expected. We conclude that the cutting prescription for the two-loop web stated in 
eq. (4.48) produces the correct imaginary part. The graphical representation of the cutting 
prescription (4.48) is similar to that in figure 8 , and we omit it here. 


In the above we have computed the imaginary part of the web with time-like kinemat¬ 
ics. As the web depends on two distinct angles, we may also consider the diagram in the 
case of mixed time- and space-like kinematics, for example cosh 712 = v± ■ u 2 > 0 and 
cosh 723 = — V 2 ■ V 3 > 0. !) A natural question is then whether also in this case the imaginary 
part is computed correctly by the formalism of section 3. As we shall see shortly, the for¬ 
malism readily applies, with the one difference that the imaginary part has no contribution 
from the V-cha 11116 !, as propagators stretched between mutually space-like eikonal lines 
have vanishing cuts. Put differently, the discontinuity in the V’-channel does not contribute 
to the imaginary part, cf. eq. (3.14). 

Returning to the formula for the imaginary part in the explicit form (4.49), we observe 
that in the above case of mixed time- and space-like kinematics, the roots of the propagator 
P 23 (y) lie outside the range of integration, cf. the remarks below eq. (3.16). As a result, 
the second term in eq. (4.49) vanishes and in the first term the principal-value prescription 
may be dropped, 

c ^ 

Im-^i 2 1 ) = / dx d y A i2 (x) Pjj ] (y) log - . (4.59) 

v ’ ’ Jo y 

9 The opposite-type kinematics Vi ■ V 2 < 0 and V 2 ■ V 3 > 0 is of course equivalent by the antisymmetry of 
the web under the interchange 712 < —► 723 of the cusp angles. 
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After inserting eq. (4.52) and integrating out the delta functions, the imaginary part is 
readily expressed in terms of multiple polylogarithms; these in turn can be simplified into 
classical polylogarithms, yielding 

Im = Ji2(x)W)(-2G( / o 2 (^),0; 1) + 2G(p 1 tyO,0; 1) 

+ G(0; p 2 {x))G(p 2 ('tp); 1) + G(0;pi(x))G(p 2 (YO; !) 

- G(0;p 2 {x))G(pi('tp); 1) - G(0; pi(x))G(pi(ip); 1)) 

= 7r-R(x)-R(^)(Li2(^) -logif>(\]ogil> -log(l - VO +log(x + 1) - glogx) - C 2 ) • 

(4.60) 

Using eq. (4.57) it is straightforward to verify that this result agrees with the imaginary 
part acquired by eq. (4.55) upon the analytic continuation x -9 — 1/x ~~ We conclude 
that the formula (3.12) reproduces the correct imaginary part of the web also for mixed 
time- and space-like kinematics, as expected. 

5 Position-space cuts of eikonal diagrams with internal vertices 

In this section we turn to the application of the formalism of ref. [1] to diagrams with 
internal (i.e., three- and four-gluon) vertices. Here we provide details on the calculation 
of the imaginary part of the diagram involving a three-gluon vertex connected to three 
Wilson lines, as illustrated in figure 12. 



Figure 12. The non-planar two-loop three-gluon vertex diagram. 


The integrated result for the diagram in figure 12 was first obtained in refs. [34, 35] 
using a Mellin-Barnes representation of the two loop-momentum integrals. In terms of the 
cusp angles 7 ij, defined through cosh 7 ^ = —Vi ■ Vj, it is given by 

9 2 3 

F 3g = -if abc T?T|Tg - (g) Y, £i > k ^ coth ^ ‘ (S- 1 ) 

i,j,k=l 

This expression is valid for an unphysical configuration with space-like kinematics for all 
pairs of Wilson lines, i.e. Vi ■ vj < 0, as indicated by the tilde on In agreement with our 
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observations in section 2 , F*9 has no imaginary part. In contrast, in a physical configuration 
of massive Wilson lines, each velocity is constrained to the unit three-hyperboloid, either 
inside the future light cone or inside the past light cone. There are two inequivalent physical 
configurations, shown in figure 13. In both configurations at least one pair of Wilson lines 



Figure 13. Physical configurations for three distinct Wilson-line velocities: (a) one incoming and 
two outgoing lines, or (b) three outgoing lines. Configurations related by time-reversal are omitted. 
In this section we compute the imaginary part of the three-gluon vertex diagram in the physical 
configuration (b). 


is time-like separated, i.e. Vi ■ Vj > 0 , leading to a non-vanishing imaginary part. 

In the following we consider the configuration in figure 13(b), where all Wilson lines 
correspond to outgoing states, such that Vi ■ Vj >0 for each pair of Wilson lines. The 
analytic result for such time-like kinematics, denoted by F 3g , is obtained from the space¬ 
like expression F 3g in eq. (5.1) by analytic continuation 7 —> m — 7 ^ for all i / j (cf. 
eq. (3.18)). The imaginary part of the resulting expression is 

ImF 3g = -f abc T a 1 T b 2 T c 3 (g) - Y, ( 7 ^ 7 fc* - 2vr 2 7ij ) coth 7fci . (5.2) 

i,j,k =1 

In the remainder of this section our task is to compute this imaginary part with our 
formalism. As in section 3 we first need to extract the leading divergence of F 3g in its 
position-space representation. The three-gluon vertex diagram has a only a simple pole in 
e. This divergence is extracted from the radial integral over the three-gluon vertex position. 
Having extracted the leading divergence, the diagram can be written as ^ x (finite). We 
then apply position-space cuts to the finite function. The one-dimensional integrals along 
the Wilson lines are then trivially performed using the delta functions arising from the 
cut. The remaining integrations over the direction of the three-gluon vertex are performed 
numerically, after which the final result is compared to the analytic expression in eq. (5.2). 
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We start by writing down the position-space representation of It reads 10 


771 fdbcrr\CLrr\b rr\C ^ , Ae 

r S » = -/ T^Ts^J 


d D 


X 


r 4—6e 


3 / q \ 

^2 £ v k Vi ■ v i Ci C k i-T^r g((i, e)J g(Cj, e) y(Cfc, e) . 

i,j,k =1 ^ 4 ' 


(5.3) 


Here the three-gluon vertex position x is integrated over all of Minkowski space. In the 
integrand x is decomposed into a radial distance r and direction u, via x^ = ru such 
that tt 2 = 1 for time-like x and u 2 = — 1 for space-like x. Dot products between u and the 
Wilson line velocities are denoted by Ci = v i ' u - The one-dimensional integrals along the 
Wilson lines are contained in the functions g{Ci, e ), which are defined as 



dxj _ 

2 Xi Ci + x 2 ) + ir]} 1 ~ e 


(5.4) 


After a change of variables to hyper spherical coordinates in eq. (5.3), the radial integral 
contains the overall divergence 11 and may easily be performed, yielding a factor of 4 (. 
Restricting attention to the leading order in e allows us to set e = 0 in the finite function 
T 3 g, yielding 


fi, = -r*° TfT‘TI (. 


J~'iq — 


[~i 2 d3u £ ijk v i-VjCiCk(^rg(Ci,Q)]g{Cj,d)g(Ck,o) ■ (5.5) 

Jm ' i j h , J 


The integration domain for the three-gluon vertex direction u is MP ’ = H[j_ U Hi U dSj^ U 
dS^, the union of the upper and lower sheets of the unit three-hyperboloid and three- 
dimensional de Sitter space, defined by 


H(j_ = {«£ M 1,3 : u 2 = 1 and u$ ^ 0} , 
dS^ = {«£ M 1 ’ 3 : u 2 = — l and uq ^ 0} . 


Having written the three-gluon vertex diagram in the form in eq. (5.5), we are ready to 
apply our formalism to obtain the imaginary part of T 3g from its cuts. This in turn gives 
the imaginary part of the full diagram F 3g , as they are proportional up to a real constant. 


The imaginary part of F 3g is computed from the formula in eq. (3.12), as illustrated 

schematically in figure 14. The cut propagators stretching between the three-gluon vertex 

and the Wilson lines take the obvious form Ai(x*) = —tt5(u 2 — 2 XiCi + x 2 ), rather than 

eq. (3.10) for propagators connecting two Wilson lines. In order to resolve the support 

- 1,2 

of these delta functions in the different subregions of MP it is convenient to introduce 
variables y* that are equal to Ci, possibly up to a sign depending on the location of u. 

10 See ref. [80] for the corresponding position-space representation of this diagram in Euclidean space. 

11 This divergence is regulated by including an exponential damping factor in the radial integral (cf. 
eq. (2.3)). 
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Figure 14. Graphical representation of position-space cuts applied to a diagram with an internal 
vertex. 


Explicitly, we let m = (i for u £ U dSjj~ and y% = —(,% for u £ H 3 U dS^~. In this way Q 
flips sign between the ± regions, but m does not. In terms of these variables a cut operator 
acting on a function g(Ci,0) yields 


Cut Xi g((i, 0) = -7T / d Xi S(u 2 -2xi(i + xf) = 
Jo 


-|(i± i) (v? - 1) 


- 1/2 


- 1/2 


in H± 
in dSg , 
(5.7) 


while the principal-value part, due to Cut x . g{Ci,0) with i ^ j, evaluates to 


PV 


f 


d Xi 


_ ( ± arcosh(yj) (yf - l) 
o — (u 2 — 2 Xi (i + x\ ) + ig |=p arsinh(yj) ( yf + 1) 


- 1/2 

- 1/2 


in H^. 
in dSg . 


(5.8) 


The right-hand side of eq. (5.7) shows that the cut vanishes for u in the region Hi. This 
is a consequence of the delta function having no support inside the domain of integration 
[0, oo). In the regions dS^ and there are respectively one and two solutions to the 
delta-function constraint, as can also be understood by inspection of figure 15. Focusing on 
the contributions to the imaginary part from dS^ and dSg , we see that both the single and 
the triple cuts acting on the product g(Q, 0) g((j, 0) g((k, 0) in F^ g yield the same results in 
both regions. But apart from this product of g' s, the sum in eq. (5.5) also contains CiCk 
which differs by a sign between the two regions. As a result, the imaginary part arising 
from the regions dS^ and dSj cancel each other. The upshot is thus that the imaginary 
part of J~ 3 g arises solely from the region Hi . 


The final step in the computation of Im JFzg is now to perform the integration over u, 
the direction of the three-gluon vertex. We do not have analytic results for the integrals in¬ 
volved, but a numerical evaluation is sufficient to show agreement with the analytic formula 
in eq. (5.2). Let us give a few details regarding the setup of the numerical integration. 

The three-gluon vertex direction u may be parametrized explicitly in terms of Minkowski 
angles ijj, $ and cf>. As discussed above, the imaginary part arises solely from the region Hi, 


- 34 - 




Figure 15. Spacetime pictures of a Wilson line along the positive cc°-axis (with normalized velocity 
Vi) and a three-gluon vertex u, located in the region (a) H 1 !, (b) dS^ or (c) H+. The region dS^" is 
very similar to dSj(" and is omitted. Each figure shows two lightlike gluons emanating from u: one 
along the future light cone and the other along the past light cone. Of these two on-shell gluons 
respectively zero, one or two gluons are able to connect to the Wilson line. In other words, in case 
(a), (b) and (c) there are respectively zero, one and two solutions to the delta functions coming 
from the cut operators. This means that in the region H_ the operator Cut^ vanishes, thereby 
producing no imaginary part. Furthermore, as discussed in the main text, the imaginary parts from 
the space-like regions dSlj" and dS^" cancel each other. The only contribution to the imaginary part 
of the diagram thus arises from the region Hi. 


which may be parametrized as 



' u° = cosh^> 
u 1 = sinh ijj sin d cos </> 
u 2 = sinh ip sin '0 sin <fi 
k u 3 = sinh ip cos $ , 


0 < vp < oo 

0 < ‘d < 7 T 

0 < cf> < 2tt 


( 5 . 9 ) 


To facilitate the numerical integration over vp E [ 0 , oo) we perform a further change of 
variables z = tanh '0, which has the effect of producing a finite integration domain z E [ 0 , 1 ]. 
Explicit expressions for C* = Hi = u ■ V{ in terms of z, the angles i9, (j) and the cusp angles 
7 ij are obtained by choosing a convenient Lorentz frame. For example, 


< = ( 1 , 0, 0, 0 ) , 

= (cosh712, 0 , 0 , sinh712) , ( 5 . 10 ) 

U3 = (cosh713, 0, sin03 sinh713, cos $3 sinh713) . 


These velocities manifestly satisfy vf = 1 and v\ ■ 77 = cosh 71 j for j 
identity, ■ v% = cosh 723, fixes 03 in terms of the cusp angles, 


cos 03 


cosh 712 cosh 713 — cosh 723 
sinh 712 sinh 713 


2 , 3 . The remaining 


( 5 . 11 ) 


The explicit parametrization of the Wilson-line velocities in eq. ( 5 . 10 ) breaks the antisym¬ 
metry of J-^g under interchange of any pair of cusp angles at the integrand level. However, 
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Figure 16. (Color online.) A comparison between the numerical and analytical results for the 
imaginary part of the three-gluon vertex diagram. The left pane shows Im as a function of 
two cusp angles from numerical integration (in blue) and from the analytical result (in light gray) 
superimposed. The numerical result fits the analytic function rather well, with an overall scale 
factor deviating from one by about 2 • 10~ 5 . In the right pane we observe that A, the absolute 
difference between the numerical and the analytical values, (in red) is below the numerical errors 
(in translucent green) for nearly all points ( 7125723 ) and is on average about three times smaller. 
The relative difference with respect to the analytic formula, A/| Im |, is on average of the order 
of 2 percent, leading to the conclusion that there is excellent agreement between the numerical and 
analytical results. In both plots 713 = 0.5, while 712 and 723 vary between 0.0 and 1.2. 


the antisymmetry must be recovered after integration (cf. eq. (5.1)). At the level of nu¬ 
merical integration this indeed happens for small cusp angles, while for large cusp angles 
numerical instabilities arise from the integration near z ~ 1, i.e. very large 0- Averaging 
over the cusp angles 7 ^ remedies those instabilities. 

Having constructed expressions for 0 in terms of z, the angles $, <j) and the cusp angles 
7 ij, the imaginary part of the three-gluon vertex diagram is explicitly given by 

r l d z z 2 f 2rr 3 ( d \ 

Im^ ~ 2 , 2 / dosin'# / d</> Vf e ijk v t ■ Vj (i Ck ww- £({0}) , (5.12) 

Jo {t-z) Jo Jo i ^L 1 \ a( ,i J 

where the cut operators are absorbed into the function f?({Ci}), given by 

£({0}) = (Cut Xi + Cut Xj . + Cut Xfc - Cuta; i)X ^ )Xfc ) g(Ci, 0 ) 5 ( 0 ,0) 9{Ck, 0) 

7 r arcosh 0 arcosh Ck + ^ arcosh 0 arcosh Cfc + vr arcosh (i arcosh (j — ^ r 3 

= (Cf -1) 1 / 2 (C? -1) 1/2 (Cl -1) 1 / 2 ■ 

(5.13) 

For the numerical integration of eq. (5.12) we have used GSL [81]. A comparison between 
the numerical and analytical results for Im JAg is shown in figure 16. We find that the 
relative difference between the numerical and analytical results are at the percent level, 
with absolute differences smaller than numerical errors. We conclude that our formula for 
the imaginary part in eq. (5.12) is in excellent agreement with the analytic expression in 
eq. (5.2). This suggests the applicability of the formalism introduced in ref. [1] to obtain 
the imaginary part of any eikonal diagram with internal vertices. 
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6 Conclusions 


In this paper we have provided algorithms for the compution of the position-space cuts of 
eikonal diagrams introduced in ref. [1] and discussed the interplay of the cutting prescrip¬ 
tion with non-Abelian exponentiation. The cutting prescription is applied directly to the 
position-space representation of an eikonal diagram and computes its imaginary part to the 
leading order in the dimensional regulator e. The prescription is stated in eqs. (3.10)-(3.12). 
The relation of the imaginary part to the branch cut discontinuity is given in eq. (3.14). 

Momentum-space cuts of eikonal diagrams, analogous to the Cutkosky rules for stan¬ 
dard Feynman diagrams, were introduced in ref. [21] where they were used to show that 
the exchanges of Glauber-region gluons (i.e., maximally transverse gluons) produce imag¬ 
inary parts of the Wilson-line correlator. Any given momentum-space cut separates the 
eikonal diagram into two disjoint subdiagrams, putting the eikonal and, depending on the 
cut, possibly also a number of standard Feynman propagators on shell (see figure 5 for an 
illustration). As a result, momentum-space cuts have the conceptual advantage of factoring 
eikonal diagrams into on-shell lower-loop and tree diagrams which can be computed as inde¬ 
pendent objects. In practice, however, the resulting cut diagrams involve integrations over 
two-, three-, four-, ... particle phase space. The evaluation of these phase-space integrals 
poses a substantial computational challenge, limiting the applicability of momentum-space 
cuts for computing imaginary parts. 

In contrast, position-space cuts do not factor the eikonal diagram into disjoint subdi¬ 
agrams, but rather constrain the gauge bosons exchanged between the energetic partons 
to be lightlike. For space-like external kinematics such exchanges are causally impossible, 
and the imaginary part vanishes. For time-like kinematics such exchanges are allowed and 
generate a nontrivial evolution of the phases of the parton states, leading in turn to a close 
relation of the imaginary part of the cusp anomalous dimension to the static interquark 
potential. Position space thus offers a causality viewpoint on the origin of the imaginary 
part of the eikonal diagram. This is complementary to the unitarity viewpoint provided 
by momentum space—i.e., that the imaginary part arises from the hard partons going on 
shell and exchanging Glauber-region gluons. At the computational level, the number of 
position-space cut diagrams contributing to the imaginary part of a given eikonal diagram 
is in practice smaller than the number of momentum-space cut contributions, and several 
of the position-space cut diagrams can be seen to be equal a priori. 

We have applied our formalism to several two- and three-loop eikonal diagrams, find¬ 
ing agreement with results previously obtained in the literature [21, 28, 35, 38, 39]. These 
computations also serve to demonstrate that the position-space cut diagrams contributing 
to the imaginary part of a given eikonal diagram can be evaluated in practice in non¬ 
trivial cases. In particular, for eikonal diagrams without internal vertices— i.e., QED-like 
diagrams—the contributing cut diagrams can be evaluated systematically by means of our 
algorithm for computing the principal-value integrals involved (supplemented with a slight 
generalization of the algorithm of that in ref. [76] for expressing multiple poly logarithms 
in terms of ones with constant indices). 

The formalism developed in this paper allows us to compute the imaginary part of 
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eikonal diagrams only to the leading order in e. This appears to limit the applicability of the 
approach, but in practice Wilson line correlators can often be expressed in terms of diagrams 
with simple poles in e (once the diagrams are expressed in terms of the renormalized 
coupling). This owes to the non-Abelian exponentiation theorem [4, 5, 8-14] which states 
that the correlator can be expressed as the exponential of specific linear combinations of 
diagrams mutually related by permutations of the soft-gluon attachment points. These 
linear combinations, called webs, have the property that the leading divergence of the 
constituent diagrams cancels, leaving in many cases webs with simple poles in e. The 
organization of the exponent of the Wilson line correlator in terms of webs is particularly 
beneficial for the applicability of the present cutting prescription: the cuts must be applied 
to an entire web rather than its constituent diagrams separately, as the separate imaginary 
parts, computed to leading order in e, will cancel. In this sense the cutting prescription has 
a nontrivial interplay with non-Abelian exponentiation. 

It would be intriguing to investigate whether the position-space cuts studied in this 
paper can be utilized, or serve as inspiration, for developing (generalized) unitarity meth¬ 
ods [48-52, 54-57] for correlators of Wilson lines. Another interesting direction for future 
research is the extension of the present formalism to computations of imaginary parts of 
Wilson line correlators to subleading orders in e. 
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A Real and imaginary parts of multiple polylogarithms 

The position-space cut prescription in eq. (3.12) produces principal-value integrals when 
applied to eikonal diagrams beyond one loop. In practice, we compute such integrals as 
the corresponding full integral (which evaluates into multiple polylogarithms) minus its 
imaginary part, cf. eq. (2.10). In this appendix we describe how to construct the required 
real and imaginary part of multiple polylogarithms to arbitrary weight in a systematic way. 
Explicit formulas for imaginary parts are given up to weight four, while the real parts are 
obtained by subtracting the imaginary part from the original function. 
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Let us first introduce some notation. Multiple polylogarithms are defined recursively 
by 

f x d t 

G(ai,... ,a n ;x) = / -- G(a 2 , ■ ■ ■, a n ; t) for (ai,..., a n ) / 0 n , (A.l) 

Jo t~ a i 

starting from the special cases 

G(; 0) = 0, G(;x) = 1, G(0 n ] x) = — log n x , (A. 2 ) 

n\ 

where a n = (a,..., a) denotes a vector with n equal indices. Multiple polylogarithms satisfy 
a variety of properties. They form a shuffle algebra, 

G(ui,..., a ni 5 x) G{ a ni+l 5 • • • j a n\+n 2 ! x ) — E G^(®(r(l)) ■ ■ • j ®<r(ni+n 2 )! ■ (A.3) 

<TGE(ni,ri2) 

They are invariant under a common rescaling of all arguments: setting a = (oi,..., a n ) we 
have 


G(k a; k x) = G(a-, x) for o, / 0 and k E C* . (A.4) 

They reduce to classical polylogarithms in certain cases, 

G(0 n ; x) = 1 log n x , G(a n ; x) = — log n (l - -) , 
n! n! V a/ 

G(0„_i,a;x) = -Li.„ (^j , G(0 n , a p ; x) = (-l) p S n>p , (A.5) 

or to harmonic poly logarithms, introduced in ref. [82], 

G(a;x) = (—l) k H$(x) if Va* G a: ajG{±l,0}, (A.6) 

where fc denotes the number of +l’s in a. We refer to ref. [74] for more details. 

Having set the notation, we turn to the problem of constructing the real and imaginary 
parts of multiple polylogarithms. As it turns out, it is most convenient to obtain the real 
part as the difference of the full function and its imaginary part, 

Re G(a\ x) = G(a; x) — i Im G(a; x). (A.7) 


Thus, it suffices to determine the imaginary part. In the remainder of this section we 
accordingly focus on the construction of imaginary parts. 

The imaginary part of a multiple polylogarithm arises when one or more of its indices 
ai are located along the path of integration of the corresponding iterated integral, that is 
ai £ [0, x\. By giving either the argument or the indices an infinitesimal imaginary part, the 
imaginary part of any multiple polylogarithm is fixed recursively in terms of the imaginary 
part of the classical logarithm. This is most easily seen in the special case of a polylogarithm 
with all indices equal to zero by using its definition in eq. (A.2) in terms of logarithms, 


Im G( 0 n ; x ± irj) = — Imlog^lx ± irj) = — Im (log |x| ± mOi—x))' 
n! n\ L 


(A.8) 
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The imaginary part on the right-hand side may be obtained by simply expanding the prod¬ 
uct and collecting the terms proportional to i. The other two special cases in eq. (A. 2) are 
real constants and thus have a vanishing imaginary part. This concludes the computation 
of the imaginary part of all the special cases listed in eq. (A.2). 

We thus turn to determining the imaginary part of a multiple polylogarithm in the 
general case G(a\x) with at least one non-zero index, cf. eq. (A.l). Since an imaginary 
part arises when some of the indices a* are located along the path of integration, it will 
be necessary to know the relative locations of the a^s in the complex plane. To this end, 
we will make a few assumptions on the indices and arguments of G(a;x). First, let us 
observe that all polylogarithms encountered in this paper will have real indices a* (up to 
an infinitesimal imaginary part whose sign is fixed by the Feynman rules). In addition, 
we will assume, without loss of generality, that the last index is non-zero, and that the 
endpoint of integration x is real and positive. 

The fact that the latter two assumptions may be imposed without loss of generality 
follows from the properties of multiple polylogarithms. Indeed, multiple polylogarithms 
with any number of trailing zeros may be expressed, with the help of the shuffle algebra 
in eq. (A.3), in terms of multiple polylogarithms with a non-zero last index, multiplied by 
pure logarithms. For example, 

G(a, 0; x) = G(a ; x)G(0; x) — G(0, a ; x ), 

G(a, 0, 0; x) = G(a; x)G( 0, 0; x) — G( 0, a, 0; x) — G( 0, 0, a; x) 

= G(a; x)G( 0 , 0; x ) — G(0, a; x)G(0; x) + G( 0 , 0, a; x) . (A.9) 

Their imaginary part is thus given in terms of the (real and) imaginary parts of multiple 
poly logarithms with a non-zero last index and pure logarithms. The latter are known from 
eq. (A.8). Thus, it suffices to determine the imaginary part of polylogarithms with a non¬ 
zero last index. 

Now, taking the last index to be non-zero, we may apply the rescaling relation eq. (A.4) 
with k = —1, 


G(a; x) = G(—a; — x) for a n ^ 0, (A.10) 

to map any negative argument x to a positive argument [83]. Likewise, a complex argument 
can be mapped to a real number after rescaling by k = 1/x, yielding G(a/x ; 1). 

Our task is thus to determine the imaginary part of a multiple polylogarithm in the 
general case G(a\ x) with a non-zero last index, and with the endpoint of integration x 
being real and positive. Starting at weight one, we let a, x E M, a / 0 and x > 0. Explicit 
computation of the imaginary part—using the fact that the imaginary part of an integral 
equals the integral over the imaginary part—yields 

f x d t [ x 

Im G(a =t iry; x) = Im / - ; - = ±7r / dt 5(t — a) = ±ir 8(a)8(x — a). (A.11) 

Jo t — (adz irj) J Q 

This is the explicit formula for the imaginary part at lowest weight. 
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Moving on to higher-weight poly logarithms, let us consider a weight n multiple poly- 
logarithm. The integrand takes the form 


X = 


h 


1 1 

(ai ± irj) t -2 - (a 2 ± irj) 


1 

tn ~ ifln ± if]) ' 


(A.12) 


The imaginary part of X may be broken into real and imaginary parts of each of the 
above fractions. The imaginary part of a single fraction is simple, because it localizes the 
corresponding integration variable to a point. In contrast, the real parts do not simplify, so 
it is convenient to have as few real-part evaluations as possible. To this end, we split the 
imaginary part of the integrand (A.12) into products of imaginary parts and real parts of 
products by recursively applying 


Im(PQ) = Im P Re Q + Re P Im Q , 

Re P Re Q = Ke(PQ) + Im P Im Q , (A. 13) 


where P and Q represent either a single fraction or products of fractions. For example, 


Im(a6) = Im a Re 6 + Im 6 Re a , 

Im (abc) = Im a Re(6c) + Im 6 Re(ac) + Im c Re(a6) + 2 Im a Im b Im c, 

Im(abcd) = Im a R e(bcd) + Im b R e(acd) + Im c R e(abd) + Im d R e(abc) + 2 Im a Im b Im c Re d 
+ 2 Im a Im b Im d Re c + 2 Im a Im c Im d Re 6 + 2 Im 6 Im c Im d Re a, (A. 14) 


where a, 6 , c and d represent individual fractions in the integrand X of a multiple polylog¬ 
arithm. Notice that each term on the right-hand side of eq. (A. 14) contains only one real 
part (of a product) and a product of imaginary parts (of single factors). 

The imaginary part of higher-weight multiple polylogarithms are then computed by 
applying eq. (A. 14) to the integrand and integrating out the delta functions arising from 
eq. (2.10). For example, the imaginary part of a weight-two multiple polylogarithm for 
a, 6 , x e M, 6 / 0 and x > 0 is computed as follows, 


Im G(a ± ir], b ± ii 7 ; x) 


= dt I (hi Im 

Jo Jo 

rx rt 


1 


1 


= / dt du 

Jo Jo 

rx 


Im 


t — (adz ir]) u— (6 ± irj) 

1 \ _ / 1 


t — (adz ir /) 


Re 


u — (6 dz irj ) 


+ Re 


1 


= dzn / dt 5(t — a) du Re 

Jo J 

= dzTT 0(a)6(x — a) Re 


1 


, dz 7 r / dt Re 
u-(b±irj)J J 0 

du „ f x dt 


t — (adz irj ) 
1 


Im 


t-(a±irj)J J 0 


u — (bdz irj) 
du 8(u — b) 


dz 7rRe 


— 6(b)0(t-b) 


Jo u — (b dz irj) Jo t — (adz irj) 

= dzTT 0(a)6(x — a) Re G(b ± irj ; a) dz rr 0(b)Q(x — b) Re (G(a ± irj-, x) — G(a dz irj ; 6 )) . 


(A.15) 

The real parts of lower-weight functions on the right-hand side are known inductively from 
expressions for lower-weight imaginary parts, together with eq. (A. 7) for the real part. We 
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remark that when a = 0, the first term in eq. (A.15) vanishes because G(b ± ir]]0) = 0. 
The imaginary part of weight-three and -four multiple polylogarithms have been computed 
along the same lines. Before quoting the results, it is advantageous to introduce some 
notation, in terms of which the formulas assume a nice form. 

The imaginary parts are conveniently expressed in terms of two new functions: a slightly 
different notation for multiple polylogarithms, together with generalized step functions. 
Let us first introduce the former, which is an iterated integral where the base point of 
integration may be freely chosen, 

rdn+l 

I (&05 j • • • ? ^ n '■> Q"n-\- 1 ) = / ~ I (&05 1 • • • i with I (oo5 5 1 • 

Jao t — a n 

(A.16) 

Setting the base point to zero we obviously recover the multiple polylogarithms defined 
in eq. (A.l), up to a conventional reversal of the indices. Although eq. (A.16) appears to 
define a larger class of integrals, it actually does not. Any / can be written as a linear 
combination of (products of) G' s. This is achieved by splitting the range of integration into 
a difference of paths with basepoint zero, cf. ref. [75]. For example, 

7(a 0 ; ai; a 2 ) = 1(0] oi;a 2 ) - 1(0; ar, a 0 ) = G(ai\a 2 ) - G(ar,a 0 ) ■ (A.17) 


At higher weight one splits the innermost integrals first. At weight two, 
pos d t 

i ao t 02 


r a 3 dt 

7(a 0 ; ai, a 2 ; a 3 ) = / -- (G(ai]t) - G(ai; a 0 )) 

J an t~a 2 K 


= G(a 2 ,oi;o 3 ) - G(a 2 ,ai;a 0 ) - G(ai; a 0 )(G(a 2 ; a 3 ) - G(a 2 ;a 0 )) . (A.18) 


In this way any I can be written in terms of C s. 

We also introduce generalized step functions 9(a \,..., o n ), which may be thought of 
as enforcing a\ < ■ ■ ■ < a n . In terms of ordinary single-variable step functions, 

n— 1 

6(cLi , . . . , 0-77,) — n 9(a,i -|_i — at) for n > 1 and Vi: a* 7 ^ aj+i and aj G 1. (A.19) 
2=1 

Equal adjacent arguments are dealt with using the following definition 

6(.. ., a, b,.. ., b, c,... ) = — -9(. .., a, b, c,...) for real indices a / b / c. (A.20) 

n times 

Infinitesimal imaginary parts produce an overall sign, 

9(..., a ± iri,...) = ±9(..., a,...) for atR. (A.21) 
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In this notation, the imaginary part of multiple polylogarithms up to weight four are 


Im G(a\x) = t r 0(0, a, x), 

Im G(a, b;x) = tt 9(0, a, x) Re G(b] a) + tt 0(0, b, x) Re /(6; a; x ), 

Im G(a, b, c; x) = tt 9(0, a, x) Re G(b, c; a) + n 9(0, b, x) Re [ I(b ; a-x)G(c-b)] 
+ tt 9(0, c, x) Re /(c; b, a; x) + 27t 3 0(O, c, b, a, x), 

Im G(a, b, c, d; x) = tt 9(0, a, x) Re G(b, c, d\ a) 

+ tt 9(0, b, x) Re [ I(b ; a; x)G(c, d; &)] 

+ tt 0(0, c, x) Re [ I(c ; b, a; x)G(d ; c)] 

+ tt 9(0, d, x) Re I(d; c, b, a; x) 

+ 27t 3 0(O, c, b, a, x) Re G(d ; c) 

+ 27t 3 0(O, d, b, a, x) Re I(d; c ; b) 

+ 27t 3 0(O, d, c, a, x) Re /(c; 6; a) 

+ 27t 3 0(O, d, c, b, x) Re I(b; a; x), 


(A.22) 
(A.23) 

(A.24) 


(A.25) 


where a, b, c,d,x 6 M, x is positive, and in each case the last index is non-zero. The 
suppressed Feynman irj ’s may be reinstated by replacing the indices according to a —> a±ir/. 
We conclude this section by providing proofs of eqs. (A.23) and (A.24). 


Proof of eq. (A. 23). We proceed by direct computation, using eq. (A. 14) and the identity 

Im JBo = =F7 t<5(0> 


ImG(o, b ; x) = Im 


nx rt 

I dt / du 
/ t—a—iO I u—b—iO 

'0 Jo 

rt 


= / d t dulm 

Jo Jo 

rx rt 


1 


1 


t—a—iO u—b—iO 


= dt du 

Jo Jo 

rx rt 

= dt du 

Jo Jo 


Im 


(mo) Re + Re (^) Im (^io) 


ir5(t - a) Re ( + Re ( u ~ b ) 


= tt Jo dt § ( t - a ) J o duRe + n J o *Re(ffio) J o du s(u-b) 

= vr0(O, a, x) f duRe(^ I zF=lo) +7r / dt Re (i=^io) 0 ( 0 ’ 6 ’ *) 


= vr0(O, a, x) Re 


du 

u—b—iO 


+ 7T0(O, b, x) Re 


dt 


t—a—iO 


= tt9(0, a, x) Re G(b] a) + 7T0(O, b, x) Re I(b] a; x). 
This completes the proof of eq. (A.23). 


(A.26) 
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Proof of eq. (A. 2f). We proceed by direct computation, using eq. (A. 14) and the identity 
Im ?Eo = z F 7r5 (0> 

Im G(a, b, c; x) = Im OTo J q ^=io / s4=io) 

f*t 


r*x nt pi 
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Jo Jo Jo 


du Im f 


i i 
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\ t—a—iO u—b—iO v—c—iO 


= dt du du 

Jo Jo Jo 
+ Im 
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t—a—iO 


Re 
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u—b—iO v—c—iO 
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+2 Im 


u—b—iO 
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t—a—iO 


Re 


Im 


l _ l 

t—a—iO v—c—iO 


+ Im 
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v—c—iO 


Re 


l _ l 

t—a—iO u—b—iO 


1 

u—b—iO 


Im 


l 


v—c—iO 


f*x rt (* u 

= Jo dt Jo dU Jo dV “ a) Re ( u-ft-ffltt-c-io ) + * 6 ( U ~ 6) Re ( t-a-iOv-l-iO 

+7T<5(u - c) Re + 2n 3 5(t - a)8(u - b)8(v - c) 

= ^( 3 , 1 ) + 1(3,2) + -^( 3 , 3 ) + 1(3,4) ■ (A.27) 

We continue by evaluating each term in the last line of eq. (A.27) separately. The first term 
evaluates to, 


t*X /*t f*U 

I(3,i) = 7T dt S(t — a) du duRe 
Jo Jo Jo y 

= 7:6(0, a, x) j du f duRe ^ 1 ,.^,^ ) 


/o JO 

r*a t*u 


= 779(0, a, x) Re du J du^^^^ 
= 770(0, a, x) Re G(b, c; a). 


The second term in eq. (A.27) evaluates to, 

t*X f*t f*U 

/ (3,2)=7r/ dt d u8(u-b)l dvRefj^a^ 

Jo Jo Jo v 


(A.28) 


du Re 


t—a—iO v—c—iO 


( dt 9(0, b, t ) 

Jo 

9(0,b,x) dtj^ duR 


—a—iO v—c—iO 

1 1 


= n0(0, b, x) Re dtj^_ ^ du^^Q 
= 770(0, b, x) Re (i(b; a; x)G(c\ b )) . 


(A.29) 
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The third term in eq. (A.27) evaluates to, 


rx ft ru 


1(3,3) = * J d t J d u J dv 5{v - c ) Re ( t _^_ i0 u J b _ i0 ) 
= 7T Io dt L du 6 , (°’ c ’ n ) Re (rooffio) 

= 7T lo dt6 ^' C ^ f c ( t—a—iO u—b—iO ) 

= tt0(O ,c,x) d tj' duRe( t _^_ i0u _l_ i0 ) 

= 7r0(O, c, x) Re d du t _'_ i0u _l_ i0 ^ 

= 7T0(O, c, x) Re I(c; b , a; x). 


The fourth term in eq. (A.27) evaluates to 

pX ft fU 

zn = 


px pt pu 

( 3 , 4 ) = 27t 3 / dt du d v 8(t — a)5(u — b)5(v — c) 

Jo Jo Jo 

px pt pu 

= 27r 3 / d t 5{t — a) / du5(u — b) / dv 5{v — c) 
Jo Jo Jo 

= 27r 3 [ d t 5{t — a) f du 5(u — b)6(0,c,u) 

Jo Jo 

= 2ir 3 / dt — a)9(0, c, 6, t) 


(A.30) 


J o 

= 27r 3 0(O, c, b, a, x). (A.31) 

Adding up the four contributions in eqs. (A.28)-(A.31) according to eq. (A.27) we obtain 
the result given in eq. (A.24). This completes the proof. □ 

The identity (A.25) may be shown by completely analogous steps, and we therefore omit 
its proof here. 


B Algorithm for achieving canonical-form polylogarithms 

In their position-space representation, eikonal diagrams without internal vertices, and the 
corresponding cut diagrams, take the form of iterated integrals. A first step in the evaluation 
of these diagrams therefore amounts to recognizing the definition of multiple polylogarithms 
in terms of iterated integrals, cf. eq. (A.l). The resulting multiple polylogarithms depend 
in general on the kinematical variables of the problem through both their indices and their 
argument. They may be rewritten in terms of polylogarithms with constant indices, in so- 
called canonical form, by exploiting the Hopf algebra structure of multiple polylogarithms, 
which encodes the plethora of functional relations within this class of functions. In this 
appendix we describe the algorithm to cast multiple polylogarithms in canonical form, 
which is extensively used in the computations in section 4. 
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The algorithm to cast multiple polylogarithms in canonical form relies on the Hopf 
algebra structure of multiple polylogarithms. In particular, we make use of the notions 
of the symbol and coproduct of multiple poly logarithms, see refs. [74, 75] and references 
therein, as well as a procedure from ref. [76] to map symbols to poly logarithms. In order to 
describe the algorithm, we start by setting up notation. Let us denote an eikonal diagram, 
or a cut eikonal diagram (or a partially integrated result thereof), by the function g(x), 
depending on n variables, x = {x \,..., x n }. (This set of variables typically contains the cusp 
angles and possibly some remaining integration variables.) After recognizing the definition 
of multiple polylogarithms, the function g{x) is given in terms of functions G(a(x); z(x)) 
depending on x through both their indices a(x) and their argument z(x). To simplify the 
following presentation, we will assume 12 that g(x) has uniform transcendental weight w. 

The algorithm follows three steps. Let us first state the algorithm and subsequently 
elaborate on each of the steps separately. 

Algorithm (Multiple polylogarithms in canonical form) 

1 . Compute the symbol 5[gf(T)] of the function g(x). 

2. Apply a map Mg to the symbol <S[<7(a?)], whose purpose is to construct a polyloga- 
rithrn in canonical form with the same symbol as the original function. The resulting 
expression differs only from the original function by terms proportional to transcen¬ 
dental constants (which are in the kernel of the symbol map). 

3. Compute subsequently the coproducts A 2 ,i i ... ) i, ..., A„,_ii to reconstruct 

any missing terms proportional to constants with transcendental weight 2, 3,.. ., w— 1 
respectively. 

The output of this algorithm is a new function h(x) in canonical form, which is numerically 
equal to the original function g(x). In the remainder of this appendix we shall give the 
definitions of the symbol S, the map Mg and the coproduct A. We conclude by illustrating 
the application of this algorithm to the non-planar two-loop ladder diagram considered in 
section 4.1. 

The first step of the algorithm involves computing the symbol of multiple polyloga¬ 
rithms, first introduced in ref. [84]. The idea of the symbol map is to encode the functional 
relations among multiple polylogarithms as simple algebraic identities in the corresponding 
tensor algebra. 

One way to define the symbol is by considering the total differential, 

n-1 / _ \ 

dG(a n -i ,..., cii, o n ) — ^ ' G(a n _i, ■ ■ •. Oj_|_i, flj—i, ■ ■ ■, oi, fl«) dlog [ - - ) , (B.l) 

“ \Oj Oj_l J 

12 If this assumption fails, the algorithm may be applied to each subexpressions of uniform weight. 
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and to define the symbol of a multiple polylogarithm analogously, cf. ref. [72], 


n— 1 

S\G((l n — i, . . . , &n)] = ^ ^ S , . . . , Q-i-|-l, CLi— 1 , • • • 5 &l5 &n)] 

Z— 1 


/ 0-2 & Z +1 

\&z &z-l 


(B.2) 


in the case of generic indices a*; i.e., non-zero and mutually different. The formula for the 
symbol in eq. (B.2), augmented with formulas for special cases and the rules of symbol 
calculus (for which we refer the reader to refs. [74, 75]), allows the symbol 5[g(x)] to be 
computed. This completes the first step of the algorithm. 

The second step of the algorithm takes the resulting symbol as input and maps it to an 
expression of multiple polylogarithms in canonical form, whose symbol is the same as the 
symbol of the original function. A procedure that achieves this goal was given in appendix 
D of ref. [76]. We cast their procedure in the form of an explicit map and furthermore make 
a slight generalization in order to deal with functions of more than two variables which 
have a sufficiently factorized form. Let us first define the map and then point out wherein 
the slight generalization resides. 

The map M$ is defined recursively in the number of variables. Starting with the case 
of a single variable, we define the map M x which takes tensors to functions, 


M X (T) = X ) liT= ( a ™ x + h ^ ® •" 0 {aix + bi) , ^ 

I 0 otherwise , 


where a* and bi are independent of x. The map M x is linear in the space of tensors: 
given a symbol S = ^TcjTj, with rational numbers c- L and tensors T t . one has M X {S) = 
])T Cj A'l x (T t ). The map M X (T) is designed to construct a function in canonical form, such 
that its symbol is given by T plus possibly other tensors which have at least one entry 
independent of x. The proof of this statement was given in ref. [76]. 

Generalizing to the multivariate case, we let x = {xi,..., x n } denote a collection of at 
least two variables, and define the multivariate map 



where the projectors P map functions to functions according to 

p£l(h) = h + M£\S-S[h}) . (B.5) 

The maps M^\s) which occur on the right-hand sides of eqs. (B.4)-(B.5) are linear in the 
space of tensors (as in the case of a single variable) and are defined to act on elementary 
tensors T by recursion in the number of variables, 




' M X1 (T) 

< M X2 ^ Xn (T^ ® • • • ® T«) M X1 (T ( r+1 ) ® • • • ® TW) 
k M X2 ,..., Xn (T) 


for r = 0 , 
for r = 1, 2,.. 
for r = w . 


w — 1 , 


(B.6) 
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For the map to be non-vanishing, the last w — r indices must depend on x\. 

Its output is then given by a canonical function with argument x\ and weight w — r, 
multiplied by an xi-independent function of weight r. The projectors P^'g(h) add such 
functions to their input, thus gradually constructing a function in canonical form, starting 
with functions of x\ with weight w down to weight 1 and repeating the process for the 
remaining variables X 2 , X 3 , ..., x n . As a result, the multivariate map M$(S) generates a 
function of the form 

Y c h,-,iu G (ai n ;x n ) ■ ■-Giai^Xi), (B.7) 

where the 5i k are independent of xi,... , x *,■ This expression is by definition in canonical 
form. For a single variable, the indices are in fact constants and as slight abuse of terminol¬ 
ogy this is what we occasionally refer to as canonical form, bearing in mind that eq. (B.7) 
is the proper definition of canonical form in the multivariate case. 

Having defined the multivariate map, let us point out the difference with respect to the 
procedure described in ref. [76]. A slight generalization resides in the definition of 
In particular, two of the maps on the right-hand side of eq. (B.6) depend on all remaining 
variables X 2 , ■ ■ ■ ,x n , rather than just the next variable x' 2 - This alteration allows us to 
reconstruct functions of more than two variables which have a sufficiently factorized form. 

A simple example illustrates the point. Consider g(x,y,z) = log x log y log z . This 
function is already in canonical form, but the algorithm should nonetheless be able to 
reconstruct this function from its symbol. In this case the weight is w = 3, and we set 
x = {x, y, z}. The symbol of g(x, y, z) is given by 

S = x®y®z + x®z®y + y®x®z + y®z®x + z®x®y + z(&y®x. (B.8) 

Let us apply the map to this symbol, 

Afe(S) = (eg O eg o Pg o Mg) (S) . (B.9) 

Since S contains no tensors with all three entries depending on x, the first map gives zero, 

Mf\s) = M X (S) = 0 . (B.10) 

Subsequently, the first projector acts on this result. Inserting its definition from eq. (B.5), 
it reduces to applying the map Ml 1 \s), which is non-zero only for tensors whose last two 
entries depend on x. Since S does not contain such tensors, the result is zero, 

P£\0) = M£\S) = 0 . (B-11) 

A non-vanishing contribution is found in the next step, coming from the two tensors which 
contain x in the last entry. 

Pf\ 0) = Mf\s) = M y , z (y ® z) M x (x) + M y , z (z ® y) M x {x) . (B.12) 

From the definition in eq. (B.3) we have M x (x) = G(0; x) = logx. The other map M y z is 
computed along the very same lines, with the results M y , z (y <g) z) = 0 and M ViZ {z ® y) = 
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( 3 ) ( 2 ) 

logy log z. By now we have that M$(S) = Plg(h), with h = Pi ( 0) = log x log y log z. 

/o\ 

Because S = S[h\, the last projector Pi ; becomes the identity, resulting in 

M$(S) = log x log y log 2 = g(x,y,z) , (B.13) 

as required, because our algorithm should reconstruct the original function from its symbol. 
This example shows that it is crucial to have both variables y and z as parameters in the 
maps My z in eq. (B.12). If one would use a single variable only, then the vanishing maps 
M y {z®y) = M y (y®z) = 0, as one can easily verify from eq. (B.3), would lead to a vanishing, 
and hence incorrect, result. This illustrates the purpose of the slight generalization of the 
map. 

In conclusion, the second step of the algorithm constructs from the symbol S[< 7 (a;)] an 
expression in canonical form, which has the same symbol as g(x), by applying the map 
M s (S) in eq. (B.4). 

The resulting expression differs from the original function g(x) only by terms propor¬ 
tional to transcendental constants, because such terms are in the kernel of the symbol map. 
Finding these missing terms is the task of the next step in the algorithm. 

The third step in the algorithm revolves around the coproduct of multiple polyloga¬ 
rithms, which generalizes the concept of the symbol. Before describing how it may be used 
to construct missing terms proportional to transcendental constants, let us first define the 
required coproducts of the form A p , q ,..., r - They are derived from the general coproduct A, 
which is defined by its action on iterated integrals with a freely specified base point of 
integration, see eq. (A. 16), cf. ref. [71] 

A [7(00, ® 1 ; • • • j Q>m ®n+l)] — 

n k 

E E I (do 5 ,..., cii k 5 a n -\- 1 ) ® IPK; a ip+ 1 ) ■ ■ • 1 a i P +l — li a i P + 1 ) • 

k =0 0=io<i\<-~<ik<ik+i=n+l p =0 

(B.14) 

The right-hand side of eq. (B.14) consists of tensors with two entries, each entry having a 
weight between 0 and n, such that the two weights add up to n, the weight of the original 
function . 13 All possible pairs of weights are thus (0, n), (1, n — 1),..., (n, 0). Grouping the 
tensors by those pairs of weights decomposes the coproduct into 

A = X] A m . (B.15) 

p+q=n 

In other words, the action of A Pi9 on /(ao;ai,... ,a n ;a n+ 1 ) yields the subset of terms in 
eq. (B.14) of tensors with weight (p,q). Besides the operator A Pi9 , the third step of our 
algorithm also uses operators with multiple indices A pq r . Those are defined recursively 
in terms of A p q . For example, A M>r is defined by application of A qjr to the second entry 

1,i The weight of I(ao;ai,... ,a„;a„+i) is equal to its number of indices, n. Likewise, a pair of weights 
(p, q ) is attributed to a tensor T p (g) T q where the weights of T p and T q are p and q , respectively. 
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of all tensors obtained from A Ptq+r . These are the definitions of the coproducts which are 
needed in the third step of the algorithm. For more details, see ref. [75]. 

The coproducts thus defined may be employed to construct missing terms proportional 
to transcendental constants, starting with the lowest-weight constants, cf. ref. [76]. To be 
specific, in the previous step of the algorithm we constructed a function h(x) such that 
S[g(x) — h(x)] = 0. This means that the difference g(x) — h(x) must be proportional to 
transcendental constants. Following the recipe in the third step of the algorithm we act 

on this difference with the coproduct A2.1.1 to first find terms proportional to (2- The 

coproduct takes the form 

A 2 ,i„..,i [g(x)-h(x)\ = ^ (x) ® log R l?j (x) (g) ■ ■ ■ <g> log R lw (x), (B.16) 

where Li 3 ^ is a linear combination of weight-two multiple poly logarithms, and the 
Ri 3 ,..., Ri w are rational functions. Since the weight-two object Li 3i ___^ w must be propor¬ 
tional to £2, we write Li 3 ^ = kC ,2 for some rational number k. This constant of pro¬ 
portionality can be determined by numerical evaluation at some specific values for x using 
Ginac [85] and running the PSLQ algorithm [86, 87]. The hereby obtained transcendental 
constant k C2 multiplies a polylogarithmic function, whose symbol is given by Ri 3 ®- ■ ■<S)R , lvj , 
arising from the tail of the arguments of eq. (B.16). Feeding this symbol back into the first 
step of this algorithm and collecting the resulting multiple polylogarithms from the output 
of step two produces a function h 2 (x) in canonical form, which is to multiply the constant 
2- As a consequence we have that 

A 2 ,i,...,i [g(x) - h(x) -k(- 2 h 2 {x)) = 0, (B.17) 

and we conclude that the difference g(x) — h(x) — k (2 h 2 (x) is equal to terms proportional 
to transcendental constants of weight three and higher, which are in the kernel of Ay 1. 
Iterating this procedure with the coproducts Asy ...y, ..., A^-yi allows us to moreover 
reconstruct the other missing constants with transcendental weight 3,— 1, respec¬ 
tively. The final output of the third step in the algorithm is thus a rewritten version of the 
original function in canonical form. 

Let us conclude this appendix by demonstrating an explicit application of the al¬ 
gorithm, involving in particular the use of the coproduct. To this end, we consider the 
expression in the first line of eq. (4.19), J = 2R(x) 2 g(x)> where the interesting part is 
given by the weight-three function 

ff(x) = ^ (G(p i, 0, p\\ 1) — G(pi,0, p 2 \ 1) - G(p 2 , 0, pi; 1) + G(p 2 , 0, p 2 ; 1)) . (B.18) 

It depends on a single variable x through its indices, p 1 = and p 2 = jz~z- We wish 
to express this in terms of multiple polylogarithms with constant indices and argument %. 
Following the algorithm, we start by computing the symbol of g(x), 

S[g(x)} = X ® (1 - x) ® X - X®X®X + X ® (1 + x) ® X • (B.!9) 
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Application of the map from step two yields a multiple polylogarithm in canonical form 
with the same symbol, 

h(x) = M x (S[g(x)]) =G( 0 ,-l, 0 ;x)-G( 0 , 0 , 0 ;x) + G( 0 ,l, 0 ; X ) . (B. 20 ) 

Indeed, the symbol of the difference vanishes, S(g(x) — h(x)) = 0. Yet the functions g(x) 
and h(x) are not equal, because they differ by terms proportional to zeta values up to 
weight three, in this case £ 2 and £ 3 - According to the third step of the algorithm, the 
missing terms proportional to £2 are obtained first by acting with the coproduct A 2.1 on 
the difference, 


a 2 ,i [g{x) ~ h{x)] = c 2 <g> G( 0 ; x) = A 2 ,i [c 2 G{ 0 ; x)} , (B. 21 ) 

where we expect that c 2 = fc £ 2 for some rational number k. Explicitly, we find 

C 2 = 5 ^; 1) + \G{^ 1) - Li 2 (x) - Li 2 (- X ) + |Li 2 (»ji) 

- |Li 2 (l - x) ~ log(x) log(l + x) + \ log 2 ( X ) - log(l - x) log(x) ■ (B.22) 

Evaluating this expression at any value of x with Ginac yields 

c 2 = -0.822467033424113218... = -±£ 2 . (B.23) 

Inserting this result for c 2 into eq. (B.21), we conclude that A 2X \g{x)~ (h{x)—^(2G(0] \'))] = 
0. To find the last missing contribution proportional to £ 3 , it suffices to evaluate the re¬ 
maining difference numerically at a single point, 

g(x) - (Mx) - ^G( 0;x)) = -0.601028451579797142... = -^£3 . (B.24) 

We finally conclude that 

g(x) = Hx) - |£ 2 G'(°; x) - |£3 

= G( 0 , - 1 , 0 ; X ) - G( 0 , 0 , 0 ; x ) + G(0,1,0; x) - |£ 2 G(0; x ) - ^£3 • (B.25) 

We have thus succeeded in expressing g(x) in terms of multiple poly logarithms with con¬ 
stant indices and argument x- Inserting this result into T 1 ' 2 ' 1 = 2R(x) 2 g(x) reproduces 
the second line of eq. (4.19). It is now a simple matter of applying eq. (A.5) for multiple 
polylogarithms with constant indices, to rewrite g(x) i n terms of classical polylogarithms, 
thus reducing the expression to the form given in the last line of eq. (4.19). 

This completes the illustration of our algorithm in a practical example and thereby also 
completes our description of each of the three steps in the algorithm to rewrite multiple 
polylogarithms in canonical form. 
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